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ABSTRACT 

We examine the disc-jet connection in stellar mass and supermassive black holes by in- 
vestigating the properties of their compact emission in the X-ray and radio bands. We compile 
a sample of ^100 active galactic nuclei with measured masses, 5 GHz core emission, and 2- 
10 keV luminosities, together with 8 galactic black holes with a total of ~ 50 simultaneous 
observations in the radio and X-ray bands. Using this sample, we study the correlations be- 
tween the radio (Xr) and the X-ray (Lx) luminosity and the black hole mass (M). We find 
that the radio luminosity is correlated with both M and Lx, at a highly significant level. In 
particular, we show that the sources define a "fundamental plane" in the three-dimensional 
(logL R) logLx, logM) space, given by log Lr = (0.60±°;^) log£ x + (0.78l°;^) logM + 
7.33l 4 Q7, with a substantial scatter of or = 0.88. We compare our results to the theoretical 
relations between radio flux, black hole mass, and accretion rate derived by Heinz & Sunyaev 
(2003). Such relations depend only on the assumed accretion model and on the observed radio 
spectral index. Therefore, we are able to show that the X-ray emission from black holes accret- 
ing at less than a few per cent of the Eddington rate is unlikely to be produced by radiatively 
efficient accretion, and is marginally consistent with optically thin synchrotron emission from 
the jet. On the other hand, models for radiatively inefficient accretion flows seem to agree well 
with the data. 

Key words: accretion, accretion disks - black hole physics - X-rays: binaries - galaxies: 
active - radio continuum: general - X-rays: general 



1 INTRODUCTION 

The ultimate observational evidence of a celestial body being a 
black hole comes from dynamical studies, by measuring the grav- 
itational influence of the central object on neighboring stars and 
gas, both in binary systems and in the nuclei of galaxies. However, 
there are a number of distinctive signatures of black hole-powered 
activity that are usually regarded as proxy of black hole existence. 
Relativistic jets emitting synchrotron radiation in the radio band are 
one such signature, the second most common being the presence 
of strong, compact power-law X-ray emission commonly associ- 
ated with the inner part of an accretion flow 1 . Indeed, as a general 
property, accretion onto compact objects and the launch of rela- 
tivistic outflows/jets seem to be correlated (or symbiotic, Falcke & 
Biermann, 1995) phenomena iBegelman, Blandf ord, & Rees 1984| 
|Rawlings & Saunders 1991} . This implies that, at some level, a cor- 
relation between jet and disc flux is unavoidable. 

Observationally, jet morphologies and spectral properties of 
both radio and X-ray cores are remarkably similar in the case of 
black holes of stellar mass (galactic black holes, hereafter GBH) 



1 For the case of stellar mass objects, these signatures are common to most 
neutron stars, too. In this cases, only dynamical studies can provide a proof 
of the black hole nature of the source. 



and of their supermassive counterparts in the nuclei of galaxies 
(hereafter SMBH). If jets are launched in the innermost parts of 
the accretion flows, as commonly assumed, then these similarities 
suggest that it should be possible to understand the physics of both 
black hole accretion and jet production by studying all those sys- 
tems as a single class. To this end, radiation emerging at radio and 
(hard) X-ray frequencies is the most direct probe of the immediate 
vicinity of the black hole: effects of extinction are unimportant in 
the radio band, while in the X-rays, where a large fraction of the 
accretion energy is radiated, such effects can be accounted for with 
good enough spectral capabilities. Thus, radio and X-ray observa- 
tions, together with the information on the central black hole mass, 
should allow us to study the relationship between the fundamen- 
tal parameters characterizing black hole activity such as the central 
black hole mass, Mbh, and the accretion rate, M, as well as the 
disc-jet coupling. 

Until less than ten years ago, quasars and active galactic nu- 
clei, and the supermassive black holes believed to power them, 
were regarded as exceptional (and extreme) objects. However, be- 
ginning with the work of Kormendy & Richstone (1995) and, in 
particular, Magorrian et al. (1998), the idea that SMBH reside in 
the nuclei of virtually every galaxy in the nearby universe has be- 
come almost commonplace. A sharp picture of the demograph- 
ics of local SMBH comes from the work of Ho, Filippenko, & 
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Sargent (1997b), who have carried out a detailed optical spectro- 
scopic survey of a large number of nearby galaxies. From this work, 
it was found that between a third and a half of the sample has 
AGN-like spectra, albeit of low luminosity, (LLAGN; either of the 
LINER, Seyfert or Transition type; see Ho, Filippenko, & Sargent 
1997a), thus confirming that SMBH are not only present in galax- 
ies, but that they are also active (at least as some level). As a mat- 
ter of fact, follow-up radio surveys of optically selected LLAGN 
<Ho & Ulvestad 20011 |Nagar et al. 2002) have yielded extremely 
high detection rates, with the radio emission having predominantly 
a compact core morphology, occasionally accompanied by jet-like 
features. X-ray studies at arcsecond (or sub-arcsecond) resolu- 
tion with the Chandra X-ray Observatory have also allowed us to 
firmly pin down the emission properties of the nearby dim galac- 
tic nuclei through systematic surveys of LLAGN l H o et al. 20011 
ITerashima et al. 2002 Ter ashima & Wilson 2003 l. These pieces of 
information on local, low-luminosity objects nicely complement 
those on more luminous (and more distant) "classical" AGN (i.e. 
Seyfert galaxies, QSOs, radio galaxies) and allow us to investigate 
the dependences of observable properties on black hole mass and 
accretion rate. 

As a matter of fact, alongside the work on black hole 
demographics, there have been various claims for the existence 
of a correlation between radio luminosity or radio loudness 
(a measure of the ratio of radio to bolometric luminosity) 
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has also led to the suggestion that there is a systematic difference 
in radio loudness between the two classes, SMBH being on 
average more radio loud iFalcke & Biermann 1996 1. Furthermore, 
among X-ray binaries, black holes tend to be more radio loud than 
neutron stars (Fender & Kuulkers 2001 1, possibly also indicating 
a mass dependence of the radio loudness parameter. However, 
the difficulty of separating the dependence of the radio power 
output on the accretion rate (due to the lack, or the neglect, of an 
independent observational indicator of it for SMBH) unavoidably 
makes the evidence of any such correlation rather weak. Indeed, 
other authors, using different samples, have recently found no 
evidence for such correlations <Ho 2002al |Woo & Urry 2002") 
|Oshlack, WeEster , & Whitt ing 2002) . 

A link between the disc accretion rate and the generation of 
relativistic radio jets has been suggested by Willott et al. (1999) 
for radio galaxies and by Ho & Peng (2001) for Seyfert 1 nu- 
clei, on the basis of the strong observed correlations between ra- 
dio and optical powers. Moreover, by studying a sample of galactic 
nuclei with measured black hole masses, Ho (2002) found clear 
evidence of radio loudness being anti-correlated with some estima- 
tor of the dimensionless accretion rate. Interestingly, Brinkmann 
et al. (2000) have also correlated the ROSAT All-Sky Survey and 
the VLA FIRST catalog, and have shown that there is some degree 
of correlation between the monochromatic X-ray luminosity at 2 
keV and the radio (5 GHz) luminosity for bright AGN and Quasars 
(for both radio loud and radio quiet sources), while Sambruna et 
al. (1999) have found a weak correlation between lobe radio power 
and 2-10 keV luminosity in a sample of radio loud AGN. 

The situation is rather different for GBH where the mea- 
sured black hole masses span a very narrow range (around 10 
solar masses). In this case, the observed, rather large, changes 
in luminosity can only be explained if they are somehow linked 
to variations in the accretion rate. On this note, recent work has 
shown that GBHs in the low/hard state (where their spectra are 



dominated by a power-law X-ray emission with typical bolomet- 
ric luminosities Lboi/iEdd ^ 10 -2 ) always possess compact radio 
cores (Fender 2001 1 the radio luminosity of which is tightly cor- 
related with their X-ray luminosity over more than three orders 
of magnitude, with Lb oc L ^ 7 (Ga llo, Fender, & Poole y"2"002| 
|Gallo, Fender, & Pooley 2003| ICorbel et al. 2003> . The evidence 
of such a correlation emphasizes the effects of the dependence of 
the radio luminosity on the accretion rate and hence the connec- 
tion between accretion and jet processes. However, an alternative 
model has been put forward recently, based on the broad band 
spectral energy distributions of some GBH in the low/hard state, 
in which the X-ray emission from those objects is produced by 
optically thin synchrotron radiation emerging from the jet itself 
(Markoff, Falcke, & Fender 2001 1. The tightness of the observed 
correlation between radio and X-ray luminosity, and its slope, seem 
to support such an idea ( Markof f et al. 2003IICorbel et al. 2003L 

Finally, we should also note that GBH with bolometric lumi- 
nosity close to their Eddington limit display a more complicated 
behavior when observed in the radio band. There is evidence that 
in the high/soft state (when the spectral energy distribution is domi- 
nated by a quasi-thermal component with kT ~ 1 keV) continuous 
jet production is inhibited, while at still higher luminosities (in the 
so-called very high state) powerful, episodic super-luminal ejection 
events have been observed (in particular in the prototypical micro- 
quasar GRS 1915+105, Mirabel & Rodriguez, 1994). Thus, high 
luminosity black hole X-ray binaries seem to display a (temporal) 
dichotomy between radio loud and radio quiet states, reminiscent 
of that of powerful quasars. 

This rapid (albeit somewhat disorderly) progress in the obser- 
vational studies of black holes of all classes seems to offer the pos- 
sibility for understanding fundamental scalings of black hole prop- 
erties with Mbh and M. Accretion disc theory can provide us with 
the relationships between the observed radiative output (hence X- 
ray luminosity) and the black hole masses and accretion rates. Un- 
fortunately, theoretical modeling of the relation between jet proper- 
ties (related to the observed radio luminosity) and the physical pa- 
rameters of the system (Mbh, M) has been hampered by the lack 
of a clear understanding of the mechanism by which jets are accel- 
erated and collimated. Relating the X-ray and radio properties of 
a black hole system to only Mbh and/or M has not been straight- 
forward. However, in a recent paper, Heinz & Sunyaev (2003) have 
demonstrated that, under the general assumption that the jet forma- 
tion process is not qualitatively different among SMBH of different 
mass or between SMBH and GBH, it is in fact possible to derive a 
universal scaling between the jet (radio) luminosity at a given fre- 
quency, L v , and both mass and accretion rate. The derived relation 
is independent of the jet model and has scaling indices that depend 
only on the (observable) spectral slope of the synchrotron emission 
in the radio band, and on the accretion model (see §|5|below). 

Motivated by these findings, the aim of this work is twofold. 
First, we want to examine (or re-examine) the significance of the 
correlations between radio and X-ray luminosities with Mbh and 
with each other in a large sample of black holes with measured 
masses that includes both GBH and SMBH. This will be dealt with 
in the first part of the paper: in ij2| we describe the selection criteria 
and the properties of our sample, while in Sj3|and £|4]we present the 
results of the correlation analysis and we describe the observational 
properties of the fundamental plane of black hole activity. Then, in 
the second part of the paper (J5j we will make use of the model by 
Heinz & Sunyaev (2003) to relate observed correlation coefficients 
that define the fundamental plane to theoretically predicted ones. 
This will allow us to place constraints on the physical properties of 
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the accretion flow from which the jet is launched, and in particular 
on its radiative efficiency. Section|6|is devoted to a discussion of our 
findings and of their general implications for our understanding of 
black hole activity in different regimes. Finally, we summarize our 
results in 



2 THE SAMPLE 

2.1 Supermassive black holes in galactic nuclei 

We have selected from the existing literature a sample of black 
hole-powered systems with measured masses, the nuclei of which 
have been observed both at 5 GHz (mostly with arcsecond resolu- 
tion with the VLA) and in the 2-10 keV band. For future reference 
we define the dimensionless black hole mass M = Mbh/Mq and 
accretion rate m = (Lbol/??)/iEdd = Mc? / Lsaa oc M/M, 
where r\ is the accretion efficiency. 

We first considered the full sample of ~ 40 nearby inactive, 
or weakly active galaxies with existing nuclear black hole mass 
measurements from observations of spatially resolved kinematics. 
To these we have added a comparable number of bright AGNs (and 
QSOs) with nuclear black hole mass measured from reverberation 
mapping of their broad line region (a similar sample was compiled 
by Ho et al. 2002). From this sample we selected all objects which 
have been observed in both the radio and X-ray bands. 

In order to obtain a more statistically representative sam- 
ple, we also searched the existing literature for both nearby low- 
luminosity galactic nuclei (Ho et al., 2001; Terashima et al., 2002; 
Nagar et al., 2002; Terashima & Wilson, 2003) and for relatively 
bright Seyfert nuclei (either type 1, type 2 or Narrow Line Seyfert 
1 ) and radio galaxies with available radio and X-ray flux measure- 
ments. We assign black hole masses to these systems using the 
observed correlation between black hole masses and stellar veloc- 
ity dispersion <Gebhardt et al. 20001 |Ferrarese & Merritt 207)01 2 , 
mostly using the values of the velocity dispersion provided by the 
HYPERLEDA catalogue jSimien & Prugniel 2002} . We note here 
that it is beyond the scope of this paper to discuss the different 
methods (and their qualities) for estimating black hole masses in 
the nuclei of galaxies, and refer the reader to the comprehensive 
literature on the subject (see e.g. De Zeeuw, 2003, and references 
therein). 

We do not include in our sample distant quasars, for which the 
M — a relation cannot be used because of the lack velocity dis- 
persion measurements. In practice, we want to avoid any indirect 
method for black hole mass estimation that relies on e.g. the con- 
tinuum optical/UV luminosity (as done by Woo & Urry, 2002, for 
example). 

Relativistically beamed sources (i.e. those whose jet axis 
points towards our line of sight) are dominated by the boosted jet 
emission, and cannot be used to test the disc-jet coupling. We there- 
fore excluded from our sample BL Lac objects. Among the Quasars 
in our sample, only 3C 273, which has an extremely high radio 
loudness and a blazar-like spectrum, is likely to suffer from strong 
Doppler boosting of the radio jet. On the other hand, according to 
the unification scheme, Seyfert 2 nuclei should not be preferentially 
viewed pole on. For all the other sources (mainly low-luminosity 

2 We use here the relationship M = 1.3 X 10 8 (o-/200km s -1 ) 4 ' 58 , as 
derived in Ferrarese (2002). See however Tremaine et al. (2002) for a thor- 
ough statistical discussion of the different scalings that have been claimed 
for the M — a relation. 



AGN and Seyfert 1), for which the nature of the (relatively faint) 
radio emission is not well established, we have assumed that the 
orientation of their jets with respect to line of sight is randomly 
distributed. A more detailed discussion of the possible incidence of 
relativistic beaming as a selection effect will be presented in i|6.2l 

2.2 Galactic black hole sources 

The Galactic X-ray binaries included in our sample have been se- 
lected to have (a) simultaneous X-ray and radio observations, or 
RXTE All-Sky-Monitor (ASM) X-ray data in conjunction with ra- 
dio fluxes available from the literature, and (b) publicly available 
RXTE-ASM X-ray and Green-Bank Interferometer (GBI) radio 
lightcurves (from which we estimated the 5 GHz fluxes by inter- 
polating between the 2.25 GHz and the 8.3 GHz channels). 

We treated the GBI and ASM data of each object in the fol- 
lowing way: we first re-gridded the ASM X-ray lightcurve to the 
radio lightcurve, and then sorted the data into bins of ascending 
X-ray flux. This procedure assumes that individual bins are not 
correlated and can thus be re-ordered and rebinned to improve 
statistics. Rebinning the data in time intervals does not change 
the results significantly, but reduces the dynamic range in X-ray 
flux. In this sense, X-ray-flux-ordered binning produces a better 
representation of the scatter in this variable. Since the GBI be- 
comes noise dominated below about 15mJy, we conservatively 
chose to consider data at or below this value as upper limits and 
split the sample into detection and upper limits before binning. 
ASM fluxes were converted to 2 — lOkeV luminosities using 
the a conversion factor of 3.2 x 10 -10 [ergs -1 cm _2 ]/[cts s -1 ] 
I Grimm, Gilfanov, & Sunya ev" 2002) , which assumes a crab-like 
source spectrum. 

Black hole masses for GBH are derived from stellar kinemat- 
ics and we use the value provided in the literature. For the micro- 
quasar LS 5039, no mass estimate is available, and we assumed 
M BH = 10 M Q . 

2.3 Global properties of the sample 

Our final sample (Table0 comprises 149 entries, of which 99 are 
individual SMBH systems. The remaining 50 correspond to the ob- 
servations of 8 different GBH at various luminosity levels. Out of 
these 8 sources, 5 are so-called X-ray transients (Cyg X-3, GRO 
J1655-40, GRS 1915+105, XTEJ1 118+480 and XTE J 1859+226) 
and 3 persistent sources (Cyg X-l, GX 339-4 and LS 5039). The 
SMBH sample includes 14 Quasars; 19 Seyfert galaxies of type 
1 (in this class we include all the Seyfert spectrally classified as 
types 1 to 1.9); 32 Seyfert 2 galaxies; 7 Narrow Line Seyfert 1 
galaxies (NS1); 11 Low-Ionization Nuclear Emission line Regions 
(LINERS) of type 1.9; 13 LINERS of type 2 and 2 LINER/HII 
Transition nuclei. Separating mass estimators into direct ones (stel- 
lar kinematics, maser emission, gas kinematics and, to a lesser ex- 
tent, reverberation mapping) and indirect ones (all those that infer 
a measure of the central black hole mass from the observed Al — a 
relation), we can assess possible biases introduced in our sample 
by the different mass measurement techniques. 55 masses are mea- 
sured indirectly (4 Syl, 4 NS1, 30 Sy2, 15 LINERs and the two 
transition objects) and 44 directly (all the 14 QSOs, 15 Syl, 2 Sy2, 
3 NS1, 8 LINERs plus Sgr A* and M 32). 

For the radio and X-ray luminosities taken from the literature, 
we have assumed Ho — 75 km s -1 Mpc" 1 , correcting the quoted 
values when necessary. When available, we have included informa- 
tion on the radio spectral properties of the sources. In Table □ wc 
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Table 1. Radio and X-ray properties of black holes with measured masses 
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14 
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16 
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u 
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67 
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s 


68 
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F 
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24 
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16 
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25 
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24 
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I CT 


16 


NGC 1068 
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s 
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15 
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M 


4 
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16 
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26 
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24 
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I CT 


16 


NGC 1667 
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S2 
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40.05 


15 
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I CT 


16 
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S2 


38.99 


F 
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16 
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S 


27 
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16 


NGC 2787 
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F 


4,29 
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24,28 
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G 


10 


NGC 2841 


12.0 


L2 
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F 


29 
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28 
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I CT 


30 


NGC 2992 
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S2 
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U 


14 
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14,15 
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16 


NGC 3031 (M81) 
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4 
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F 


27,50 


40.20 


15 
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26 
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26 
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34 
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S 
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5 
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R 


21 


NGC 3362 
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S2 
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S 
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14 
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16 


NGC 3516 
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SI 
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s 
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R 


4 
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S2 


35.45 


F 


29 
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28 
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30 


NGC 3675 
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T2 
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u 


29 
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30 
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S 
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24 
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S 


4 
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R 


21 
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S2 
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s 


14 
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45 
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16 
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34 
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R 


21 
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30 
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S1.9 
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F 


4,29 


40.52 


35 
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M 


4 


NGC 4261 (3C 270) 
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L2 
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U 


11 


41.17 


24 


8.72 


G 


4 


NGC 4278 
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F 
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39.96 


31 
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S 


60 
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T2 
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U 


29 
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28 


6.80 


la 


30 
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L2 
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24 
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G 


4 


NGC 4388 


16.8 


S2 
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S 
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15 
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26 
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36 
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4 


NGC 4450 


16.8 


L1.9 


36.53 


F 


29 


40.34 


24 


7.30 


la 


26 


NGC 4457 


17.4 


L2 


<35.70 


U 


29 


39.97 


24 


6.86 


la 


26 


NGC 4472 


18.8 


S2 


36.69 


s 


27,29 


<38.80 


37 


8.80 


la 


38 


NGC 4486 (M87) 


16.1 


L2 


39.78 


F 


4,29 


40.55 


39 


9.48 


G 


4 


NGC 4494 


9.7 


L2 


<35.65 


U 


29 


38.86 


28 


7.65 


la 


30 



A fundamental plane of black hole activity 



5 



Table 1 - continued 



Tsjnp asoi 

1NLJY 4JU1 


1 A 8 
lO.o 


^9 
oz 


^A 90 
JO.ZU 


o 


97 
Z / 


AO 98 
4U.ZO 


1 A 
14 


7 00 

/ .yu 


lex 


9A 
ZO 


wrir ASA8 
INOY 4J4o 


1 A 8 
lO.o 


T 9 

IjZ 


jO.ju 


r 


90 

Ly 


30 70 
JV. ly 


31 
J 1 


7 AO 
/ .4U 


1(7 


9 A 
ZO 


Mr,r ASAS 
1MJL 4JOJ 


Q 7 


1 Q 


JO. 1 J 


r 


90 31 
Ly,3 1 


30 SA 

jy. jo 


31 
J 1 


7 70 
/ . /U 


A CT 


OA 
ZO 


wrir AS7Q 

IN UL 4 J / 7 


1 A 8 
lO.o 


1 q 


J 1 -OJ 


r 


90 31 
Ly,3 1 


A1 1A 
41.14 


98 31 
Zo, J 1 


7 SS 


Act 


9A 
ZO 


wr,r ASQA /A/1 1 Ozl'i 
1MJL 4jy4 L1V11U4^ 


Q Q 

y.o 


T 9 

LZ 


^7 9.A 


r 


A S3 

4, J J 


AO 70 
4U. /U 


1 S 


OA 

y .U4 


c 

o 


A 
4 


Mr,r AA3A 


17 
1 / -U 


T 1 Q 

L,l .y 


^A AO 


TT 
U 


90 

Ly 


/io A3 
<_ J0.4J 


37 
J / 


7 00 

/ .yu 


Ict 


38 
Jo 


wnr A79S 

1> ^ J\- 4 / Z J 


12 4 


S9 


on 


TT 

U 


Ly 


30 1A 
3y . lO 


98 
Zo 


7 AO 


Ict 


30 
JU 


wrir A73A 


A 3 


T 9 

LZ 


1A SO 


p 

r 


90 

Ly 


30 A9 

jy.oz 


9A 
Z4 


7 30 
/ . JU 


Ict 


9 A 
ZO 


wrir S033 


1 8 7 
1 o. / 


<31 S 


^A 70 
JO. ly 


c 
o 


97 
Z / 


A1 00 
4 1 . UU 


9A 98 
Z4,Zo 


7 30 
/ . JU 


Ict 


30 
JU 


wnr siqa 

1> vJV_- J 1 ;?4 


7 7 


S9 
oz 


35 50 


§ 


27 


39 80 


15 


a on 


1(7 


1 A 


"NinP S9S9 


09 3 

7Z. J 


<;9 

OZ 


^S QA 


o 


1 A AO 

14, 4y 


A9 07 

4z.y / 


1 A 
14 


8 1 9 

O.IZ 


1(7 


1 A 
lO 


Wnr 1 S973 


91 3 


«^9 


^A 99 
JO.ZZ 


p 


1 A 9A 
1 4.Z4 


^"A9 
<,4Z.U 


1 A 
14 


A 37 
O.J / 


1(T 


1 A 
1 o 


Nfif S3A7 


31 2 


S9 
OZ 


^7 i n 

3 1 . 1 U 


p 


14 55 


40 0^- 


45 


6 70 


A CT 


1 A 

lO 


Mr,r SSA8 


70 9 
/U.Z 


<31 S 


^S SS 
Jo. Jo 


c 

o 


A 97 
4,Z / 


A3 SS 
4J. JJ 


s 


8 03 
o.UJ 


p 


91 
Zl 


wrir SQ9Q 
rNOY jyzy 


33 9 


^9 
oZ 


^S ^0 
Jo. JU 


c 
o 


1 A AO 

14, 4y 


AO it 
4U. / T 


AS 
4J 


7 99 
/ .ZZ 


1(7 


I A 
lO 


Mr,r A1 AA 


1 1 A 

110 


^9 
OZ 


jy.yj 


c 
o 


AO 
4U 


AO SA 
4U. JO 


A1 
41 


Q 1 

y. iy 


1(7 


9 A 
ZO 


wr,r A9S 1 

1M UL OZ J 1 


OA 8 
y4.o 


*v9 

OZ 


AO Q1 


c 

o 


A S3 
4,J J 


A9 1 S 
4Z. 1 J 


1 s 

1 J 


8 73 
o. / J 


r 


A 
4 


wrir asoo 

iNUL OJUU 


AO 
4U.U 


T 9 

LZ 


^S QO 

jo.yu 


p 
r 


31 SA 
J 1 ,J4 


AO 1 1 
4U. 1 1 


3 1 
J 1 


8 98 
o.Zo 


1(T 


30 
JU 


wr,r 7AAQ 

IN LJ Y / 40y 


AA A 
oo.o 




^S ^8 
Jo. Jo 


c 
o 


A AO 
4,4V 


A3 31 
4 J.J 1 


A9 
4Z 


A 81 
O.o 1 


P 


91 
Z 1 


wnr 7A77 
IN / D / Z 


S3 9 


Q9 
OZ 


17 7S 

J / .ZJ 


c 

o 


1 A SS 

14, 33 


/II 17 
4 J.J / 


1 A 
1 4 


A en 

O.oU 


1(T 


1 A 
10 


wr,r 7A89 
INOY /OoZ 


AS 
DO.U 


^9 
OZ 


^S S8 
Jo.oo 


c 
O 


1 A AO 
14, 4y 


^"A3 9 
<^4J.Z 


1 A 
14 


7 9S 
/ .ZJ 


1(T 


1 A 
lO 


\Y:r 77/11 

INLjL, / /4J 


94 4 
Z4.4 


Q9 
OZ 


1A QQ 


c 
o 


97 
Z / 


1Q 71 

jy. / 1 


7A 
Z4 


A A7 
0.4 / 


1(T 


1 A 
10 


Pri 009 A-i- 1 9Q 
JrLJ UUZO+lZy 


A97 
DZ / 


r\ 
V 


AO 1 7 
f U. 1 / 


TT 

LJ 


yl 
4 


AA AA 
44.44 


A3 
4J 


7 73 


P 


91 
Zl 


pr: OOS9_i_9S 1 
JrLJ UUJZ+ZJ 1 


AGO 

oyu 




1Q A9 


T T 

U 


A 
4 


AA AA 
44. OO 


AA 
44 


e 3A 

5. J4 


p 
K 


91 
Z 1 


pr, 080A-I-7A1 
JtLj UoU'+t /Ol 


A30 
4JU 


V 


QQ A\ 


TT 
U 


/I 
4 


AA 9A 
44. ZO 


A3 
4J 


8 98 
o.Zo 


p 


91 
Z 1 


pn 08AA.I-3AQ 
x^Lj Uo44+J4y 


9A8 
zoo 




IS 1 9 

Jo.lz 


T T 

u 


/I 
4 


A3 90 

4J.zy 


43 
4J 


7 3A 
/. J4 


p 
K 


91 
Z 1 


T>Ct oqs3j.ai a 

X^Lj uyjJ-t-414 


1118 
11 lo 


n 

w_ 


AO 1A 
f U. 14 


TT 

LJ 


/I 
4 


A A SO 
44. JU 


A3 
4J 


8 9A 
o.ZO 


p 


9 1 
Z 1 


pr, 191 1 i 1 A'X 
rKJ 1Z11+14J 


3A9 
JOZ 


n 

V 


A1 OS 
4 1 .Uo 


TT 

LJ 


yl 
4 


A3 A1 
4J.0 1 


A3 
4J 


7 A1 
/ -Ol 


p 

XV 


91 
Z 1 


pr, 199A-1-093 /'3r'973~l 
rKJ IZZOtUZJ \3k_L/3) 


70S 

/UJ 


V 


AA 0^ 


p 

r 


A SA 
4, JO 


AS 70 
4J. /U 


AA 
44 


8 7A 
o. /4 


p 


9 1 
Z 1 


pn 1 99Qj_90A 
x^vj izzy-rzu4 


9A8 
ZOO 


C\ 

w_ 


^S AS 
Jo.M-J 


TT 

LJ 


/I 
4 


AA 90 

44. zy 


AA 
44 


7 88 


p 


9 1 
Z 1 


pr, 1 107-lOSS 


AQO 

oyu 


r\ 

V 


^S QS 

jo.yo 


TT 

u 


A 
4 


AA S1 
44. J 1 


AA 
44 


8 AA 
o.44 


p 


9 1 

Z 1 


pr, 1/111 i,l /i 9 
rvj 1411+44Z 


ISO 
JoU 




18 71 
Jo. / 1 


T T 

LJ 


/I 
4 


A3 S3 
4J. J J 


A3 
4j 


7 QO 

/ .yu 


p 
K 


9 1 
Z 1 


pr 1 A9Aj-01 S 
X^O 14ZOtU1J 


3AA 

JOO 


V 


^8 Q8 


TT 
U 


/I 
4 


A3 80 
4J.oy 


AA 
44 


8 A7 
O.O / 


p 


91 
Zl 


j>n 1 ai ij-ass 

rKJ IOIj+Ojo 


SAS 

JOj 




1Q 7S 

jy. / j 


p 
r 


A SA 
4,J0 


AA 3A 
44. J4 


AA 
44 


8 38 
O.Jo 


p 


9 1 
Z 1 


pr 1 700j-S 1 8 
r kj 1 /UUtj 1 o 


1 AOA 
14UO 


V 


AO Q9 


c 
o 


A S3 
4, JJ 


^-A3 90 
<^4J.ZU 


A3 
4J 


7 78 
/ . / o 


p 


9 1 
Zl 


pri 9 1 ^0-loqq 

r \j Zl JU+Uyy 


9SS 

ZJJ 


n 

V 


^8 80 
Jo.oV 


c 
o 


A SA 
4, JO 


A3 SS 
4J. JJ 


AA 
44 


8 1 A 
o. ID 


p 


91 
Z 1 


Qr 1 9n 
JY1ZU 


1 38 

1 JO 


*vi 

O 1 


A1 SS 

*+ 1 .33 


p 

Jr 


1 1 S3 

1 1,J J 


A3 QS 
4J.yj 


1 1 
1 1 


7 3A 
/ . JO 


p 


9 1 
Z 1 


Qr QQfl 1 

jy jyu.J 


9A1 
Z4-1 


O 1 


A1 00 
4 1 \)y 


c 

o 


1 1 

1 1 


AA 00 
44. UU 


1 1 
1 1 


S S3 
o. JJ 


p 


9 1 
Zl 


T ir,r A 1 00 

UUL OIUU 


1 1 £> 
110 


oz 


^8 SO 
Jo. JU 


TT 

u 


1 A 
14 


^"A3 A 
<,4J.O 


1 A 
14 


7 79 
1 . 1 L 


T 

1(T 


1 A 
1 o 


ogr A 


n ooc 




17 SO 
JZ.JU 


p 

JT 


A7 
OZ 


11 1/1 
JJ. J4 


A3 
OJ 


A A 1 
0.41 


c 
o 


A 4 
04 


r. 7(T v 1 

L,yg A-i 


009 1 


v 1 1 j n 


<^9Q AA 


P 
r 


A1 
O 1 


3A AA 
JO. 44 


A1 
O 1 


1 OA 
1 .UO 


c 
O 


AO 

oy 


r . ,., v 1 
L,yg A-l 


O Of^9 1 


/ ; i > i_i 

V J L J II 


^90 AS 

<^zy.4j 


p 
r 


A 1 
01 


3A S7 
JO. J / 


A 1 
01 


1 OA 
1.U0 


c 
O 


AO 

oy 


r, 7(T v 1 
L,yg a- i 


009 1 


v 1 1 j n 


90 

zy.oj 


p 

JT 


A1 
1 


3A AS 
J0.40 


A 1 
1 


1 OA 
1 .UO 


c 
O 


AO 

oy 


r,.„ y 1 
L-yg A- 1 


O 009 1 
U.UUzl 


/ ; i > lj 

V I L J II 


90 AS 

zy.oj 


p 

r 


A 1 
01 


3A S7 
JO. J / 


A1 
01 


1 OA 
1.U0 


c 
O 


AO 

oy 


cyg a-i 


009 1 


v j 1 3 n 


90 AA 

zy.oo 


p 
r 


A1 
01 


3A AA 
J0.04 


A1 
Ol 


1 OA 
1 .UO 


c 
o 


AO 

oy 


rVfr y 1 
L,yg a-i 


009 1 
U.UUZ1 


v 1 1 3 n 


90 9A 


p 
r 


79 
/ Z 


3A SS 
JO. Jo 


A1 
Ol 


1 OA 
1 .UO 


c 
o 


AO 

oy 


r, m v 1 
cyg a- i 


009 1 


v i [3 n 


90 30 
Ly. JU 


p 

Jr 


79 

1 L 


3A AS 
JO.OJ 


A 1 
Ol 


1 OA 
1 .UO 


c 
O 


AO 

oy 


r . ,., v 1 
cyg A-l 


O 007 1 
U.UUZ1 


/ M > II 
V I L J II 


90 1 8 

zy. i o 


p 

JT 


79 
/ Z 


3A AA 
JO.OO 


A1 
01 


1 OA 
1.U0 


c 
O 


AO 

oy 


cyg a-j 


00Q 
U.UUV 


v j 1 3 n 


^in aa 

<^ JU.04 


p 

r 


A 1 
01 


3A OA 
JO.y4 


A1 
Ol 


1 OO 
1 .UU 


c 

o 


70 

/u 


Cyg X-3 


0.009 


GBH 


31.62 


F 


61 


37.19 


61 


1.00 


s 


70 


CygX-3 


0.009 


GBH 


31.85 


F 


61 


37.41 


61 


1.00 


s 


70 


Cyg X-3 


0.009 


GBH 


32.17 


F 


61 


37.60 


61 


1.00 


s 


70 


Cyg X-3 


0.009 


GBH 


32.43 


U 


61 


37.88 


61 


1.0 


s 


70 


GR0 J1655-40 


0.0032 


GBH 


<29.76 


F 


61 


36.64 


61 


0.85 


s 


65 


GRO J1655-40 


0.0032 


GBH 


<29.69 


F 


61 


37.55 


61 


0.85 


s 


65 


GRO J1655-40 


0.0032 


GBH 


29.94 


F 


61 


35.27 


61 


0.85 


s 


65 


GRS 1915+105 


0.0125 


GBH 


<30.89 


F 


61 


38.17 


61 


1.20 


s 


65 


GRS 1915+105 


0.0125 


GBH 


<30.89 


F 


61 


38.47 


61 


1.20 


s 


65 


GRS 1915+105 


0.0125 


GBH 


<30.89 


F 


61 


38.73 


61 


1.20 


s 


66 


GRS 1915+105 


0.0125 


GBH 


31.60 


F 


61 


38.25 


61 


1.20 


s 


66 


GRS 1915+105 


0.0125 


GBH 


31.76 


F 


61 


38.40 


61 


1.20 


s 


66 


GRS 1915+105 


0.0125 


GBH 


31.76 


F 


61 


38.65 


61 


1.20 


s 


66 



6 Merloni, Heinz & Di Matteo 



Table 1 - continued 



rzv ain A 
V i A joy-** 


ft OOA 


noy 
udh 


9Q Ql 
zy .y 1 


r 


7 


36 AS 


7 


1 00 


c 






\ i A Jjy-4 


n. nn/i 

U.UU4 


uDn 


9Q G7 
Zy.o 1 


p 
r 


7 


A9 
JO.^-Z 


7 

1 


1 00 

1 .uu 


c 







r T y Tin a 


o OOA 

U.UU^t 


nou 
uDn 


90 HQ 
Zy.oy 


p 


7 


Q ^ AO 


7 


1 00 


c 


Q 



ny ian a 

V I A Jj^-f 


o ooa 

U.UU^t 


nou 
UDn 


9Q f\l 
Zy .OZ 


p 


7 


"3./^ 1 9 
JO. 1Z 


7 


1 00 

1 .uu 


c 





v i A jjy-^f 


OOA 


noy 
uDn 


9Q 71 
Ay. 1 1 


p 


7 


1 9 
JO. 1Z 


7 


1 00 

1 .uu 


Q 






rzv iin a 


ft OOA 


noy 
uDn 


9Q f\f\ 
zy .DO 


p 


7 


3fi 1 A 
JO. I'-t 


7 


1 00 

1 .uu 


c 






v i A jjy-*+ 


n ooa 

U.UUH 


noy 
uDn 


9Q A^ 
zy .^K? 


p 


7 


Jj.o 1 


7 


1 00 

1 .uu 


c 


Q 



ny Tin yi 


OOA 


noy 
uDn 


9Q 1 1 

zy . 1 1 


p 


7 


Q ^ ^9 
JJ. JZ 


7 


1 00 

1 .uu 


c 


Q 
O 


ny Qin a 


ft OOA 


noy 
UDn 


98 ^A 
Zo. J4- 


p 


7 


3A 9 1 


7 


1 00 

1 .uu 


c 



O 


r,y Tin a 
v i A jjy-^f 


OOA 
UAA/f 


noy 
uDn 


98 09 
Zo.UZ 


p 


7 


j j.y 1 


7 


1 00 

1 .uu 


Q 



O 


ny Tin /i 


n ooa 

UAAH- 


noy 
UDn 


98 


p 


7 

1 


T3 8 A 

J J.O'-f 


7 


1 00 

1 .uu 


c 


V> 
O 


ny Tin /i 


n ooa 

UAAJH 


noy 
UDn 


98 AQ 


p 


7 




7 


1 00 

1 .uu 


c 


O 


ny no a 


fi OOA 
UAA>f 


noy 
UDn 


^97 S1 
<^Z / . J 1 


p 


7 


^11 1A 
<^ J J. J4- 


7 


1 00 

1 .uu 


c 



O 


ny lin a 
v i a jjy-4- 


ft OOA 
UAAM- 


noy 
UDn 


<^Z / .ZJ 


p 


7 


j j. iy 


7 


1 00 

1 .uu 


c 



O 


LiO JUJ7 


003 


GBH 


29 80 




61 


£9 
J J.OZ 


61 


1 00 






Lo jUJy 


U.UUj 


rn it 
UDrl 


ta no 




I 


Jj.4j 


a 1 

Oi 


1 .uu 






LS 5039 


0.003 


GBH 


30.07 


s 


61 


35.67 


61 


1.00 






LS 5039 


0.003 


GBH 


30.09 


s 


61 


35.81 


61 


1.00 






LS 5039 


0.003 


GBH 


30.08 


s 


61 


35.96 


61 


1.00 






XTEJ1 118+480 


0.0018 


GBH 


<28.98 


F 


46 


35.43 


47 


1.00 


s 


48 


XTEJ1 118+480 


0.0018 


GBH 


28.92 


F 


71 


35.46 


61 


1.00 


s 


48 


XTEJ1 118+480 


0.0018 


GBH 


28.92 


F 


71 


35.57 


61 


1.00 


s 


48 


XTEJ1 118+480 


0.0018 


GBH 


28.92 


F 


71 


35.56 


61 


1.00 


s 


48 


XTEJ1 118+480 


0.0018 


GBH 


28.92 


F 


71 


35.47 


61 


1.00 


s 


48 


XTEJ1 118+480 


0.0018 


GBH 


28.92 


F 


71 


35.45 


61 


1.00 


s 


48 


XTE J1859+226 


0.011 


GBH 


<29.24 


F 


61 


36.58 


61 


>0.88 


s 


71 


XTE J1859+226 


0.011 


GBH 


<29.30 


F 


61 


36.65 


61 


>0.88 


s 


71 


XTEJ1859+226 


0.011 


GBH 


29.18 


F 


61 


36.66 


61 


>0.88 


s 


71 



NOTE: Comments: t L x calculated from the known luminosity in the 0.3-8 keV band and the observed spectral index; t Lx calculated from the GIS count 
rate, assuming T = 2 and Log(Nn) = 22.5; § L 5Ghz extrapolated from observations at 8 GHz. Col.(l): Name of the object. Col. (2): Distance in 
Megaparsecs (for H = 75 km s _1 Mpc -1 ). Col. (3): Spectral Class; GBH: Galactic black hole; L: LINER; S: Seyfert; NS1: Narrow Line Seyfert 1; T: 
Transition object (LINER/H II); Q: Quasar. Col. (4) Logarithm of nuclear luminosity at 5GHz. Col. (5): Radio spectral index (where F y = u~ aR ); F: 
flat spectrum (ag < 0.4); S: steep spectrum («r > 0.4); U: undetermined. Col. (7) Logarithm of the intrinsic rest-frame luminosity in the 2-10 keV band. 
Col. (9) Logarithm of the black hole mass. Col. (10) Mass measurement method; S: stellar kinematics; G: gas kinematics; M: maser kinematics; R: 
reverberation mapping; I[oiii] : inferred from the mass-[OIII] line-width correlation; I CT : inferred from the mass- velocity dispersion correlation. 
REFERENCES: (1) Ulvestad, Antonucci & Goodrich (1995); (2) Leighly (1999); (3) Wang & Lu (2001); (4) Ho (2002); (5) Weaver, Gelbord & Yaqoob 
(2001); (6) Nandra et al. (1997); (7) Corbel et al. (2003); (8) Hynes et al. (2003); (9) Fabbiano et al. (2003); (10) Tremaine et al. (2002); (11) Sambruna, 
Eracleous & Mushotzky (1999); (12) Pellegrini et al. (2003); (13) Saglia et al. (1993); (14) Polletta et al. (1996); (15) Bassani et al. (1999); (16) Woo & Urry 
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have marked all the sources with flat radio spectrum (qr < 0.4, 
where lvr is the radio spectral index a = —d\n L v / d\nv eval- 
uated at the frequency v =5 GHz) with flag (F). All the sources 
with qr, > 0.4 are instead classified as steep spectrum (S), while 
those for which the radio spectral index couldn't be determined are 
marked with a (U). Overall, we have 77 flat spectrum sources (32 
SMBH and 45 GBH, respectively); 43 steep spectrum sources (38 
and 5) and 29 for which the spectral slope is undetermined (all 
SMBH). 

In Figure Q we show the distributions of black hole masses, 
radio and X-ray luminosities and that of the ratio of their X-ray to 



the Eddington luminosity, Lx/iEdd = i/2-iokcV /1.3 x 10 38 M. 
It is worth emphasizing that the shape of those distributions reflects 
more the nature of our selection procedure rather than the global 
intrinsic properties of accreting black holes. For example, a precise 
measurement of the central black hole mass is a prerequisite for a 
given source to be included in the sample. As recently discussed by 
De Zeeuw (2003), not all mass ranges are equally well probed by 
the different methods, and this should introduce a strong selection 
effect in our sample. Similarly, because BH masses are more eas- 
ily measured in the nuclei of nearby galaxies, our sample is biased 
against the most luminous quasars. This should be taken into ac- 
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Figure 1. Distribution of black hole masses (top left panel), 5 GHz radio luminosity (-Lr, top right panel), 2-10 keV X-ray luminosity (Lx. bottom right 
panel) and of the ratio Lx/^Edd- Solid histograms denote detected sources, dashed ones are for upper limits. 



count when examining, for example, the Lx/LEdd distribution of 
the black holes in our sample. 

In Figure [2] we show the radio luminosity versus the black 
hole mass for objects of different spectral classes. Panel (a) shows 
the whole sample, while panel (b) concentrates on the SMBH only. 
Overplotted are some of the linear regression fits discussed in the 
recent literature (see j3|for details). 

Finally, Figure [3] shows the core radio luminosity versus the 
nuclear X-ray luminosity in the 2-10 keV band (left panel) and ver- 
sus the ratio of the X-ray nuclear luminosity to the Eddington lu- 
minosity (right panel). We represent objects in different mass bins 
with different colors to highlight a possible segregation of different 
mass bins in the Lr — Lx plane. 

In the next section we present a more quantitative statistical 
analysis of the correlations among Lr, Lx, and M. 



3 CORRELATION ANALYSIS 

Luminosity-luminosity correlations have to be tested for possible 
spurious effects introduced by their common dependence on the 
distance. Such tests can be done by performing a partial correlation 
analysis, taking distance as the third variable. Here we choose the 
so-called partial Kendall's r correlation test, proposed by Akritas 



& Seibert (1996) in the case of censored data sets. Applying this 
test to our data, we find that the Lr, — Lx correlation is strong: the 
probability of the null hypothesis (i.e. that there is no correlation) 
is less that 10~ 10 for the whole sample and ~ 1.12 x 10 -4 for the 
SMBH only (see Table|2j. 

However, it is apparent from Fig.|3|that, when the data points 
are grouped into mass bins, objects in different bins tend to lie on 
parallel tracks. Such a behavior can be seen both in the Lr — Lx 
plane (panel a), and, perhaps even more clearly, in the Lr — Lx/M 
plane (panel b). The presence of a mass segregation suggests that 
the radio luminosity of an object likely depends both on its X-ray 
luminosity and on its mass. In order to assess this hypothesis we 
once again use partial correlation analysis, now taking Lr (or Lx) 
as the dependent variable and testing its partial correlation with Lx 
(Lr) where we take M as the third variable. In addition, we also 
look for a partial correlation between Lr (Lx) with mass itself, by 
taking Lx (Lr) as the third variable. In so doing we can effectively 
discriminate between intrinsic and spurious correlations among the 
three quantities. 

The results of such tests (see Table|2j show that the radio lu- 
minosity is strongly correlated with both black hole mass and X-ray 
luminosity (P nu n < 1 x 10~ 10 for the whole data set), while in turn, 
the X-ray luminosity correlates with both mass and radio luminos- 
ity only if we include both GBH and SMBH (P nu u ~ 4.68 x 10" 5 ). 
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These results imply that any regression fit used to find correlations 
between any two variables that does not account for the depen- 
dence on the third one (as, for example trying to find the depen- 
dence of Lr on M without accounting for the dependence on Lx), 
inevitably leads to an incorrect estimate of the correlation coeffi- 
cients. 

If, in any case, we compute the correlation coefficients be- 
tween any two of these variables 3 to allow a comparison with 
earlier works, we find that logL| MBH = (29.54 ± 1.60) + 
(1.23 ± 0.20) log M, with a very large scatter (standard deviation 
4™ = 1.65). Including GBH, we find logLff = (28.75 ± 
0.18) + (1.20 ± 0.04) log M (standard deviation a\\ = 1.51), 
which is indeed very similar to what found by Nagar et al. (2002), 
but very different from what was proposed by Franceschini et al. 
(1998), based on a much smaller sample (as shown by the different 
linear regression fits plotted in Figure[2j. 

We emphasize, however, that the partial correlation analysis 
presented so far and summarized in Table[2| implies that a far bet- 
ter representation of the mutual dependencies of Lr,, Lx and M 
(and the one that minimizes the scatter) should be searched for with 
multivariate linear regression tests, that allow simultaneous fitting 
over the 3-dimensional space defined by the three variables. 



3.1 Finding the multivariate correlation coefficients 

The standard multivariate regression formalism does not treat de- 
pendent and independent variables symmetrically (see, e.g., Fasano 
& Vio 1988). In other words, if a linear regression analysis yields 
a linear regression coefficient b for the dependence of y on x, it 
does not necessarily yield a coefficient 6 _1 for the dependence of 
x on y. Closely related to this is the fact that standard \ 2 statistics, 
for which errors are only associated with the dependent variable, 
artificially reduce the estimate of the correlation coefficient in the 
presence of intrinsic scatter in the independent variable. The net 
result of these shortcomings is that x 2 distributions over the fitting 
parameter space are not fair representations of the confidence in the 
resulting best fit. 

In order to alleviate these shortcomings, we extend the statis- 
tical approach used by the "Nukers" group l Gebhard t et al. 20001 
ITremaine et al. 20 02 1 to obtain an unbiased estimator of the best 
fit regression coefficients in multivariate problems. This approach 
uses a modified chisquare estimator, called the merit function 
jFasano & Vio 19881 IPress et al. 19921 ITremaine et al. 20021 de- 
fined by 



* 2 = E- 



E,- hm. 



(1) 



where yi is the dependent variable, Xij are the independent vari- 
ables and a Vi and <r Xi . the associated error estimates, a is the zero 
intercept, and b are the linear regression coefficients, which are to 
be found. The argument of the sum in eq. Q is a measure of the 
projected distance of the data point Pi — (xij ,yi)to the regression 
hyperplane given by the equation y = a + y\ bjXj, measured in 
the x 2 space around point Pi. In other words, the argument of the 
sum in eq. Q counts the number of \ 2 intervals between Pi and 
the regression hyperplane. Because eq. Q is nonlinear in bj, it is 



3 We use a the lineal' regression method by parametric EM algorithm 
(that deals with censored data) as implemented in the ASURV package 
jlsobe, Fe igelson, & Nels on 1986} 



no longer possible to minimize \ 2 analytically. However, for a set 
of parameters b, we can still find the optimal value a m i n (b) for 
which x 2 is minimized as 



,(b) = 



(v-J2 j b i Xi 3'> 

^'<+E>^> 
Ei ( a vi +E J ( b J cr ^) 



(2) 



Finding the best fit parameters is then a simple minimization 
problem which can be treated with standard numerical minimiza- 
tion routines. Because our analysis is restricted to a bivariate prob- 
lem in this paper, it is possible to use a graphical solution in the 
form of 2-D % 2 contour plots. 

Unfortunately, such a technique cannot handle censored data. 
However, it is easy to verify that the incidence of the censored data 
points in our sample on the estimate of the linear regression coeffi- 
cient is small compared to the intrinsic scatter in the data. In order 
to do so, we have performed standard linear regression fit using 
specific techniques that are able to handle censored data, as imple- 
mented in the ASURV package jlsobe, Feigelson, & Nelson 1986) , 
on the whole data sample both with and without the upper limits. 
The differences in the resulting multivariate linear correlation coef- 
ficients in the two cases are much smaller than the estimated errors. 

In the limit of a Vi 3> <J Xij , the merit function in eq. Q re- 
covers the traditional x 2 estimator. This implies that the usefulness 
of eq. Q is limited to cases where the estimated errors are an ac- 
curate representation of the true intrinsic uncertainty, including in- 
trinsic scatter. Artificially large errors in one variable will introduce 
an asymmetry similar to that intrinsic to the traditional regression 
formalism. At face value, the errors in our sample are dominated 
by the uncertainties in the mass measurements. However, because 
we are using a sample of predominantly low redshift low lumi- 
nosity AGN, error in the distance measurements can be large due 
to peculiar velocities, which leads to large errors in the luminosity 
measurements as well. For Galactic sources, distance estimates typ- 
ically carry even larger uncertainties. Worse yet, the resulting errors 
in Lx and Lr will be correlated. As a first step, we therefore de- 
cided to follow the Nuker approach and assume isotropic uncertain- 
ties crjifj = o"l x .i = a L B ,i in all three variables and re-normalize 
these errors to produce a minimum reduced xj? cd — xV^dof of 
unity. 



4 RESULTS 

We fit the data with the function 

logL R = ^RxlogLx +£Rjv/logM + 6r, (3) 
or, alternatively with its reciprocal 

log L x = £xr log Lr + £xa/ log M + fax , (4) 

where we have used the merit function to estimate the best fit linear 
regression coefficients. Because the regression technique we adopt 
is symmetric 4 , eqs. J3J and J4j contain the same amount of infor- 
mation. 



4 We verified that the standard regression technique is instead not symmet- 
ric: for the same fitting function the best fit parameters are £rx = 0.57 
and §rm = 0.77, while by fitting eq. |4j we obtain £xr = 1-05 and 
£xm = —0.49, instead of the expected £xr = 1/£rx = 1-75 and 
£xm = — 5rm/5rx = —1.35. 
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Figure 2. Radio core luminosity at 5 GHz versus black hole mass. Upper limits are marked with arrows. Panel (a) shows the whole sample, including GBH 
(black filled circles) and SMBH, with different symbols indicating objects belonging to different spectral classes. The right panel (b) is a blow up of the SMBH 
sample. The dot-dashed line gives the regression fit proposed by Franceschini et al. (1998), the dashed line that proposed by Nagar et al. (2002), both obtained 
using different samples of SMBH only. The thick solid upper line gives the maximum core radio power as calculated by Ho (2001) for sources accreting at 
the Eddington rate. Although a correlation analysis of our sample would give results similar to those obtained by Nagar et al. (2002), as we discuss in the text, 
none of the above relations reflects the real physical scaling of radio power and black hole mass. 
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Figure 3. Radio core luminosity at 5 GHz versus X-ray luminosity in the 2-10 keV band (a) for the whole sample. Different colors correspond to different mass 
bins, while different symbols correspond to different classes of objects. In panel (b) we plot the same radio core luminosity at 5GHz vs. the ratio ^x/^Edd 
of X-ray to Eddington luminosity. The color-coding of the different mass bins makes the mass segregation more evident. 



Our final results of the multivariate regression analysis are 
shown in the form of \ 2 contour plots (using eq. |IJ as a x 2 es- 
timator) in the (£rm, £rx) parameter space in Fig.[5] The top panel 
shows the results of the multivariate fit performed on the entire 
data set. We obtain £ RX = 0.60 ± 0.11, £rm = 0.78±£'Jg, and 
&r = 7.33+to? (all the errors are one sigma confidence), with a 
dispersion <tr, = 0.88, or, cast in the form of eq. 

logL R = (0.60tE , ) ;ii)logL x + (0.78lS;;^)logA/+7.33ti^(5) 

The value we obtain for the £rx correlation coefficient is con- 
sistent, within the errors, with that found in GBH (£rx ~ 0.7) by 
Gallo et al. (2003). This also means that individual GBH sources 
for which the correlation between radio and X-ray luminosities is 
well established (GX 339-4 and V404 Cyg) do indeed follow the 



same global trend defined by black holes of all masses included in 
our sample. 

Our data set provides fairly tight constraints on the values of 
the correlation coefficients fjj — R, X,M). This in turn im- 
plies that in the 3-dimensional space (log Lb,, log Lx, log M) the 
sources are distributed preferentially on a plane, that we call here- 
after the "fundamental plane" describing black hole activity. Fig- 
ure|4|shows an edge-on view of the plane. 

The other two panels of Fig.|5]show the constraints on the two 
correlation coefficients obtained with the regression analysis on dif- 
ferent subsamples assembled according to the spectral properties in 
the radio band. For flat spectrum radio sources we obtain £rix,f — 
0.76 ± 0.13, £ RM ,f = 0.71 ± 0.10, and 6 Rjf = 1.3ll£ u £, with a 
dispersion <tr ,f = 0.81, smaller than that of the whole data set. The 
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Table 2. Results of Correlation analysis 



Variables Subsample Correlation 



X 

(1) 


Y 

(2) 




z 

(3) 


Objects 
(4) 


N 
(5) 


N x , 
ul 

(6) 


N^ 

ul 
(V) 


ul 

(8) 


r 
(9) 


a 
(10) 


Pnull 
(ID 




Log I/r 


Log 


;L X 


Log D 


All Objects 


149 


20 


14 





0.255 


0.0414 


7.3 x 10" 


10 


Log L R 


Log 


;Lx 


Log D 


SMBH only 


99 


7 


13 





0.200 


0.0518 


1.12 xl0~ 


-4 


Log L R 


Los 


:L X 


M 


All Objects 


149 


20 


14 


1 


0.448 


0.0495 


< 1 x 10~ 


-10 


Log Lr 


Los 


;L X 


M 


SMBH only 


99 


7 


13 


1 


0.450 


0.0523 


< 1 x 10~ 


-10 


Log I/r 


M 




Log L x 


All Objects 
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20 


1 


14 
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0.0469 


< 1 x 10" 


-10 


Log Lr 


M 




Log L x 


SMBH only 


99 


7 


1 


13 


0.310 


0.0547 


1.45 xlO" 


-8 


Log L x 


M 




Log L R 


All Objects 


149 


14 


1 


20 


0.184 


0.0452 


4.68 x 10 


-5 


Log L x 


M 




Log Lr 


SMBH only 


99 


13 


1 


7 


-0.022 


0.052 


0.672 





NOTE: Col. (1): Variable X. Col. (2): Variable Y. Col (3): Variable Z. Correlation between variables X and Y is studied, taking into account the mutual 
correlation of X nd Y with Z. Col. (4): subsample. Col. (5): Number of objects in the subsample. Col. (6)-(8): Number of upper limits in X, Y and Z. Col. 
(9)-(l 1): Results of partial correlation analysis, giving the partial Kendall's r correlation coefficient, the square root of the calculated variance a, and the 
associated probability P nu n for accepting the null hypothesis that there is no correlation between X and Y. 



difference in the coefficient for the Lr — Lx correlation with re- 
spect to the whole sample result may be due to a larger incidence of 
GBH in the flat spectrum subsample, which seem to have a slightly 
larger measured value of £rx jGallo, Fender, & Pooley 2002) ; the 
constraints on the other coefficient, £rm are almost as good as for 
the entire data set, and in very good agreement with it. 

On the other hand, in our steep spectrum sources subsam- 
ple, only one GBH is included (LS 5039), and this results in a 
much larger uncertainty, in particular on the radio luminosity-mass 
correlation coefficient: (£rx, s = 0.41+£^, £ RM , S = 0.94tg;?J, 
&r, s = 14.08tg gi with dispersion <t r , s = 0.87). 



5 PHYSICAL INTERPRETATION OF THE 
FUNDAMENTAL PLANE RELATION 

As stated in the introduction, a correlation between X-ray and ra- 
dio emission is expected if there is a fundamental connection be- 
tween accretion flows and jet activity. Thus, at a qualitative level, 
the existence of the fundamental plane found in Sj4|is not surpris- 
ing. However, on a quantitative level, the presence of such a plane 
and the measurement of the correlation coefficients associated with 
it provide powerful probes of accretion physics and of the inner jet. 

In the following section we will therefore lay out an avenue 
of how to use the fundamental plane relation to constrain accretion 
and jet physics. In order to do this, we shall first turn to a discus- 
sion of the theory of the radio and of the X-ray emission by the jet- 
disc system. We will first make use of the scale invariant assump- 
tion about the disc-jet coupling l Heinz & Sunya e~2003) to derive 
the expected scalings of the radio luminosity with black hole mass 
and accretion rates. Then, we will show how these theoretically 
predicted correlation coefficients translate into observable ones be- 
tween black hole mass and X-ray luminosity (£rm, £rx) for differ- 
ent models of the X-ray emission. By comparing these theoretical 
predictions to the observed values from eq. J5} we will then discuss 
the constraints we can put on the models themselves. 

5.1 Synchrotron emission from scale invariant jets 

It was recently shown by Heinz & Sunyaev (2003) that the depen- 
dence of radio luminosity Lr on black hole mass M and dimen- 
sionless accretion rate rh can be cast into a model independent form 



;/ the underlying jet physics is scale invariant. In this case, all the 
model dependent uncertainties can be absorbed into the observable 
spectral index a. The relationships between Lr and M and be- 
tween Lr and m depend only on the boundary conditions at the 
base of the jet, set by the conditions in the accretion flow feeding 
the jet. We will briefly review their argument (the interested reader 
is encouraged to consult Heinz & Sunyaev, 2003, for more details) 
before comparing their predictions to the correlations derived in the 
previous section. 

The nature and conditions in the inner disc are most likely 
governed by only a few parameters (M, rh, and possibly the black 
hole spin a). The fundamental scale imposed on the problem is 
the gravitational radius of the black hole, R e = GMbh/c 2 . let 
formation occurs in the innermost regions of accretion discs, thus 
it is natural that jet formation, too, is primarily governed by these 
three parameters only 5 . 

Then, we make the ansatz that jet structure and dynamics are 
(at least approximately) invariant under changes of R s oc M and 
m. In other words, taking a jet produced by a black hole of mass 
Mi , scaling its dimensions by a factor of M2 / Mi we should obtain 
a jet around a black hole of mass M2. 

This proposed invariance can be cast into a simple mathe- 
matical form. Any quantity / needed to calculate the synchrotron 
emission from jets can be decomposed into a structure function 
ipf(R./R s , a), which describes the spatial variation of / along the 
jet, and which depends on M only through R/i? g , and a normal- 
ization <j)f(M,rh,a), which is set by the boundary conditions at 
the base of the jet (thus, by accretion disc physics): 

/(R, M, rh, a) = <j> / (M, rh, a) ■ ^/(R/-Rg,a) (6) 

The quantities necessary for calculating the jet synchrotron 
emission are the magnetic field strength B, the jet diameter Dj e t, 
and the normalization C of the electron power-law distribution 
dtie/d^y — C^f~ p (typically, the spectral index of the particle dis- 
tribution is p ~ 2 — 3). For example, according to the prescrip- 
tion {6), the magnetic field should follow B = 4>B{M,rh, a) ■ 

5 The influence of a and the associated second scale of the problem, the 
light cylinder radius, on jet formation is unclear at this point. We will hence- 
forth keep a fixed and assume that variations in a will only introduce a 
scatter in any relation derived below which is independent of M and rh. 
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Figure 4. The edge-on view of the "fundamental plane of black hole activity". The solid line shows the best fitting function Isl . 



tpB (R/iZg, a), where <j>B — Bo = B(Ro) is the value of the field 
at the base of the jet Ro- As a geometric quantity, the jet diame- 
ter D should be directly proportional to the characteristic scale R s 
such that 4>d — Do = D(Rq) oc Ra oc M. Scale invariance also 
implies that dynamical time scales are proportional to R g /c oc M, 
and thus that characteristic velocities are scale invariant. 

Using the standard formulae for synchrotron emission 
jRybicki & Lightman 1979) , Heinz & Sunyaev (2003) showed that 
the synchrotron luminosity L v at a given frequency v emitted by 
the jet must then depend non-linearly on M and rh, following 



dln(L„) 
91n(M) 



and 

d In (m) 
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+ 
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where a is the spectral index at frequency v. 

Note that the model dependent structure functions 
ipf(R./R g ,a) scale out from these expressions. Only the 



spectral indices (a and p) and the boundary conditions 4>b and </>c 
for the magnetic field B and the electron power-law distribution 
normalization C, respectively, remain, a and p are observables: 
the electron spectral index p can be deduced from the optically 
thin synchrotron spectral index at high frequencies. The functions 
(f>B and (j>c, on the other hand, need to be provided by accretion 
disc theory. It is reasonable to assume that the relativistic particle 
pressure at the injection radius is a fixed fraction (i.e., independent 
of M and rh) of the total pressure at injection, <f>c oc cj> B . This 
leaves 4>b as the only model dependent parameter of the theory. 
Given a prescription for 4>b, we can predict how the synchrotron 
luminosity of jets should scale with M and rh. 

Since £m and £m are constants, we have in general (with 
L v = L§ GHz = Ln) 



log L R = £m log M + £m log rh + K\ , 



(9) 



where K\ is a normalization constant 6 . Variations in other param- 
eters, such as the viewing angle or the black hole spin, will only 



6 We note here that for flat spectrum jets with apt ~ 0, the canonical 
value of p ~ 2, and 4>' 2 B oc M~ x ra, the dependence of on M and 
m follows L R oc (Mrh) 17/12 = M 17 / 12 , as had been found by Fal- 
cke & Biermann (1996) for the specific case of the "canonical conical" 
JBlandford & Konigl 19791 jet model. 
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Figure 5. Shaded areas show the \ 2 density distribution, and dashed lines the \ 2 contours for the observed correlation coefficients £r,m and £r.x- The 
inner 3 contours show the formal 1,2 and 3 sigma confidence levels, the remaining contours further out show levels of Axj? od = 10. Shown are the results 
for: the entire data set (upper panel), the flat spectrum sources (lower left panel), the steep spectrum sources (lower right panel). Overplotted on each panel 
are the theoretically predicted values of the correlation coefficients where circles, diamonds and squares represent the ADAF, jet and standard disc models 
respectively. Empty symbols show the values for ajj = and filled ones for = 0.5. The lines connecting the points represent the tracks of £r,x and §r,m 
traced out by variation of or,. 
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introduce a scatter in this relationship that is independent of M and 
rh. 

Because the expressions for £m and £m do not depend on the 
shape functions tpf, they are independent of the model details. Any 
scale-invariant jet model that reproduces the observed radio spec- 
tral index cir must satisfy eqs. and J8j. This implies that (a) 
measurements of £m and £m cannot be used to constrain the func- 
tions ipf, i.e., they cannot be used to distinguish between different 
jet models but that (b) measurements of £m and £m can be used to 
place constraints on the boundary conditions at the base of the jet 
4>f, or in other words on the accretion disc model. 

While the black hole mass M is observable (see above), the 
accretion rate rh can only be inferred through radiation in other 
bands, where the emission is dominated by the accretion disc. Thus, 
in addition to a prescription of the boundary conditions (j> / through 
accretion disc theory, we must also provide a relation between the 
accretion disc luminosity (or Lx) and rh in order to model the ob- 
served I/R-m correlation. We must therefore discuss the different 
possible sources of the X-ray emission. We will start with those 
models that describe X-ray emission as produced by an accretion 
flow of some kind. 



5.2 Accretion flow origin of X-ray radiation 

For the accretion powered X-ray luminosity we can write, in anal- 
ogy to eq. |9}: 



log Lx = log M + q log rh + K2 



(10) 



where K% is a normalization constant. The efficiency coefficient q 
need not be constant, though the linear correlation analysis from ij3| 
is limited to this case. 

Using eqs. I7l9i together with eq. 1101 . we obtain the follow- 
ing general expression for the observable correlation coefficients 
as defined in eqs. J3J and 0: 



frtM 



2p + 13 + 2q r + <91n0 fl f 2p + 13 + a R p + 6a R 



p + 4 
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2p + 13 + agp + 6qr 
9 mm \ q{p + 4) 
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d In-. 



q(p + 4) 



(11) 



Different accretion models imply different values of q, and differ- 
ent scalings of the magnetic energy density at the base of the jet 
4>b with M and m, resulting in different values of d In <f>B /d In M 
and d In (f>B /d In rh. In the next two sections we calculate the ex- 
pected scalings for radiatively efficient and inefficient flows (the 
results of the calculations are summarized in Table[3}. For the case 
of radiatively efficient accretion, we will consider the X-ray emis- 
sion to be produced by a standard, geometrically thin and optically 
thick disc coupled with a hot corona. In the case of radiatively inef- 
ficient accretion flows, for the sake of clarity we will describe a de- 
tailed spectral modeling in the framework of a pure advective disk 
(ADAF). Currently this is the simplest existing dynamical model 
of a radiatively inefficient flow from which we can derive detailed 
spectral energy distributions. However, one should keep in mind 
that radiatively inefficient flows may be significantly modified by 
convection or outflows (see below). We do not, however, expect the 
arguments laid out below to be significantly modified in such cases. 

As a result of our calculations, we will show how the physical 
properties of radiatively inefficient and efficient discs do indeed im- 



ply very different scalings for the observed X-ray luminosity with 
the accretion rate rh. 



5.2.1 Radiatively efficient flows 

In the standard Shakura-Sunyaev disc model 
I Shakura & Suny aev 1973) accretion occurs via an optically- 
thick and geometrically thin disc. The effective optical depth 
in the disc is very high and photons are close to thermal equi- 
librium with electrons. The emission results in multi-color 
(T oc rh 

l/4 M -l/4 r -3/4 

) blackbody spectrum. This component 
is thought to explain the "blue bump" in AGN and the soft X-ray 
emission in GBH. However, the standard model does not predict 
the power-law X-ray emission observed in all sub-Eddington 
accreting black holes. In fact, the emission in the 2-10 keV band 
we are considering here is completely dominated by the power-law 
spectral component for SMBH and for GBH in the low/hard state. 

It is generally accepted that this hard X-ray power-law emis- 
sion is produced by inverse Compton scattering of the soft, black- 
body disc photons on a population of hot electrons that surround 
the cooler disc in the innermost region of the accretion flow (the 
so-called corona) In these models a fraction f c of the gravita- 
tional power is dissipated in the corona (and eventually emerges 
as X-ray radiation). The X-ray luminosity can then be written as 
Lx ~ fcriiLEdd oc f c rhM, and the coronal magnetic energy 
density follows B 2 oc rhf c /M (e.g.; Haardt & Maraschi 1991; 
Di Matteo, Celotti & Fabian 1999; Merloni & Fabian 2002). i x 
should therefore scale linearly with m, as long as the fraction of 
power dissipated into the corona / c does not change with the accre- 
tion rate. Note, however, that we expect some scatter in any such a 
relation due to intrinsic changes in the coronal plasma optical depth 
and temperature that may have some dependence on the accretion 
rate. 

Merloni & Fabian (2002) and Merloni (2003) have shown how 
it is possible to derive coupled equations for the accretion disc- 
corona systems under very simple assumptions about the nature of 
the turbulent magnetic viscosity inside the disc. The main property 
of their self-consistent coupled treatment is that / c is indeed con- 
stant when gas pressure dominates in the disc, so for accretion rates 
m<0.016(a v M)~ 1/8 (l - 0.84c4)" 9/8 (with a v viscosity pa- 
rameter) we should expect q = 1. On the other hand, f c oc rh' 1 / 2 
in the radiation pressure dominated parts of the disc, and q = 1/2 
(see Table[3j. It is worth noting that q and d<f>B jdrh enter eq. II II 
in such a way that the correlation coefficients £rm and £rx are 
identical for the gas and radiation pressure dominated disc cases. 



5.2.2 Radiatively inefficient flows 

The standard solution we have discussed in the previous section 
is not a unique solution for the accretion flow equations at low 
accretion rates. In recent years, much work has been devoted to 
the detailed study of low radiative efficiency accretion and in 
particular to the so called Advection Dominated Accretion Flows 
(ADAF) solutions jNarayan & Yi 1994| |Narayan & Yi 1993) 
lAbramowicz et al. 19951 |NarayanTMahadevan, & Qu^taerTT998) . 
The latter usually refers to the optically-thin ADAF branch 
which is established only for accretion rates lower than 
a critical value rh < m cr it ~ a v tReesetal. T9821 . 
However, both from the theoretical point of view 
jNarayan, Igumenschev, & Abramo wicz 20 00| and from nu- 
merical simulations jlgumenschev & Abram owicz 2000] 
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Hawley & Balbus 2002J, it has been shown that that radia- 
tively inefficient flows are prone to strong convective instabilities 
and/or powerful outflows which alter the nature of the solutions 
significantly. Despite the extensive theoretical efforts, however, the 
relative importance of convection and outflow for adiabatic flows 
is still a matter of significant debate l Balbus & Hawl ey"2002| 
|Narayan et al. 200 2 1 and no strong observational discriminant has 
yet been found. 

From our point of view, it is important to understand how the 
value of the magnetic field in such flow scales with M and rh, 
and to have a detailed model for their radiative output, in order 
to predict the scaling of the X-ray emission with these parame- 
ters. From general arguments, we expect any mechanically cooled 
flow to obey <f>% oc M _1 rh |Heinz & Sunyaev 2003} . On the other 
hand, the exact value of q can only be determined once a detailed 
assessment of the radiative processes that give rise to observed lu- 
minosity is made. This is what we discuss in the following. 

A hot accretion flow around a stellar mass (supermassive) 
black hole radiates mostly in the optical (radio) to X-ray bands. 
In the optical (radio) band, the emission results from synchrotron 
radiation. At higher energies, and up to the X-ray band, the 
emission is produced by bremsstrahlung processes for low ac- 
cretion rates and inverse Compton scattering of the soft syn- 
chrotron photons or blackbody photons from the disc in the outer 
regions when the accretion rate approaches the critical value 
{Narayan, Barrett, & M cClintock 1998 1. The predicted spectrum 
from an ADAF depends (weakly) on the ratio of the gas to mag- 
netic pressure /3, the viscosity parameter a v , and the fraction of the 
turbulent energy in the plasma which heats the electrons, 5. Here, 
we fix a v = 0.1, (3 = 10, and take S = 0.3. 

The top panel of Figure |6| shows the spectra of ADAFs (plus 
outer thin disc for rh ~ m cr i t ) for ten accretion rate values starting 
from rh = 10~ 6 up to m ~ m cr it ~ 10~ 2 . The bottom panel 
shows the relationship between the 2 — 10 keV luminosity, Lx, and 
the accretion rate, rh (for M = 10). For the chosen values of the 
parameters a v , and 8, the dependence of the 2-10 keV luminosity 
on rh is roughly given by (see also Beckert & Duschl 2002, for a 
similar calculation): 

f rh 3 ' 4 , m<10" 4 
Lx,2-iokeV oc i rh 1 ' 6 , KT 4 <rh<5 x 10" 3 (12) 
rh 2 , 5 x 10~ 3 <rh<2 x 10~ 2 

The full band spectra shown in Figure [f>| with the two vertical 
lines indicating the 2 — 10 keV energy band illustrate the origin for 
the various dependences shown in equation <12> . For rh < 10 -4 
the integrated 2 — 10 keV emission includes both the first order 
Comptonized component, which drops off at around few keV, and 
bremsstrahlung which dominates the emission at higher energies. 
At higher accretion rates the first (or higher) order Compton scat- 
tering of softer photons always dominate the 2 — 10 keV emission 
(see Figure |6| top panel). At these higher rates the optical depth 
increases and cooling becomes overall more efficient. 

This can also been seen in bottom panel of Figure |6| where 
we also plot (dashed line) the electron temperature at r — 3 as a 
function of rh. The temperature is nearly constant for m < 10 -4 ; in 
this regime the viscous heating is virtually fully balanced by advec- 
tion. At rh > 10 -4 cooling processes become more important and 
the temperature decreases. The overall change in the dependences 
in eq. 1121 reflects this behavior. 

However, it is important to keep in mind that the exact depen- 
dences in this narrow band will also depend on the microphysical 
parameters that we choose for the model. For instance, variations 
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Figure 6. Predictions from ADAF models for rh = 

io- 6 ,icr 5 ,5- 5 ,icr 4 ,5 x icr 4 ,icr 3 ,2 x icr 3 ,5 x icr 3 ,icr 2 . 

The top panel shows the broad band spectral energy distributions for a 
fixed M = 10, with the vertical dashed lines enclosing the 2 — 10 keV 
spectral energy band. For the top two models (dot-dashed lines) we also 
plot the multicolor blackbody from an outer thin disc with an inner edge at 
r = 40. In the bottom panel: the solid dots joined by the dotted line show 
the integrated 2 — 10 keV X-ray luminosity and the open symbols joined 
by dashed lines and the flow electron temperature at r = 3 versus the 
accretion rate rh. The solid line is a linear fit to the 2 — 10 keV luminosity 
over the whole of the rh with slope 2.3. 

in the electron heating parameter, S will result in changes in the 
Comptonization spectra. Therefore, in order to compare theoretical 
predictions with the observed data, the intrinsic scatter of which 
does not allow us to put tight constraints on the different emission 
regimes, we fit the Lx — rh relation over the whole range of rh 
and obtain a single power-law, Lx oc rh 2 ' 3 (i.e. q — 2.3, as shown 
by the solid line in the middle panel of Figure [6). We also calcu- 
late the relation between I/2-iokcV, versus the black hole mass, 
M for the models and, as expected, the relation is very close to 
linear with Lx oc M ' 97 . Overall the scalings derived from our 
specific computation of a specific ADAF model are consistent with 
the expected behaviour of the more general class of radiatively 
inefficient, mechanically cooled accretion flows. In fact, we note 
here that if the flow is modified by convection or outflows, we ex- 
pect the X-ray spectrum to be dominated by bremsstrahlung emis- 
sion jQuataert & Nar ayan 1999|IDi Matteo et al. 1999i . for which 
we expect q = 2. 
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Table 3. Scaling indices of the X-ray luminosity and of the magnetic field 
as functions of the physical parameters rh and M for different accretion 
flow models. 
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5.3 Jet origin of X-ray radiation 

The presence of both high energy electrons and radio/optical radi- 
ation implies that the inner jet must emit inverse Compton X-rays, 
contributing to the overall X-ray spectrum at some level. How- 
ever, under suitable conditions even the synchrotron component 
will reach X-ray energies. Especially in the presence of relativistic 
beaming, synchrotron emission could actually be responsible for 
the bulk of the X-rays. It is therefore useful to discuss the radio-X- 
ray-mass correlation expected for the synchrotron X-ray scenario 
<Markoff, Faicke , & Fend er 2001^ . 

While at low frequencies the effects of synchrotron self- 
absorption lead to the flat spectrum typically observed in core dom- 
inated sources, at high frequencies the jet becomes optically thin to 
synchrotron self-absorption, following the well know optically thin 
power-law with X-ray spectral index ax = (p — l)/2. 

Because X-rays and radio emission now originate from the 
same source, it is straightforward to derive the relation between the 
two components. The correlation coefficients £a/ and for the 
optically thick radio emission are still given by eqs. Q an d H), 
while for the optically thin X-ray emission we can simply substi- 
tute ax = (p— 1)/2 in those expressions. Finally, we can eliminate 
log rh from eq. {5) in favor of log Lx and arrive at the desired ex- 
pressions for the correlation coefficients: 

2(2p + 13 + Q R p + 6qr) 
4RX " (p + 4)(p + 5) 

z _ 2p + 13 + 2a R 6(2p + 13 + a R p + 6q r ) 

^ — vn (p+mp+5) — (13) 

where an is the radio spectral index and p is assumed to be the 
same for radio and X-ray emitting electrons 7 , and we have used a 
tilde to distinguish the synchrotron X-ray model. 

5.4 Comparison with observations 

In the previous sections we have shown how different theoretical 
models for the emission processes responsible for the observed ra- 
dio and X-ray emission from black holes can be directly translated 
into predictions for the observable correlation coefficients £rx and 
£rm (see eq. [3}. The relatively tight observational constraints on 
these indices that we have derived in section|5|can therefore be used 
to directly test these models. The theoretical coefficients for differ- 
ent models are shown in Figure [5] The circles and the squares de- 
note inefficient and efficient disc accretion respectively. Diamonds 

7 This implies that radiative cooling must either be negligible for the re- 
gion of the jet where the X-rays are produced, leaving p ~ 2 over the 
entire spectrum, or that cooling and continuous injection have produced a 
universal power-law with a slope of order p ~ 3. This limitation, while 
severe, cannot be avoided because particle transport including the effects of 
radiative cooling cannot be formulated in a simple scale invariant fashion. 



represent synchrotron X-rays from the jet. Open and closed sym- 
bols are for flat (an = 0) and steep (or = 0.5) spectrum sources, 
respectively. All model coefficients are calculated assuming p — 2 
(but see 36.41 for a discussion of steeper electron distributions as 
due, for example, to the effects of cooling). To highlight the sen- 
sitivity of the coefficients to variations in qr, we also plotted the 
tracks in £ R x — £rm space for changes in a R between and 0.5. 

For the models in which X-rays are produced by radiatively 
inefficient accretion, we have q = 2.3 and d In (f>B /d In rh = 
— d\n(f>B /din M = 0.5 (see 35.2.21 . The predicted coefficients 
for both flat and steep spectrum sources lie roughly within the 3<j 
contours of the observed sample. Thus, radiatively inefficient ac- 
cretion is consistent with the observations. This conclusion is also 
true for the general class of models for which q = 2. 

For the standard disc-corona model, we have q = 1 (high ra- 
diative efficiency) and 9 In 4>b /dlnrh = — d\n<f>B /d\riM = 0.5 
(see 35.2. 1> . The model predictions fall well outside the 3a con- 
tours of the data and are therefore inconsistent with the observa- 
tions. 

For the models in which the X-rays are produced by opti- 
cally thin synchrotron emission from the jet itself, the flat spec- 
trum model predictions are marginally consistent with the data (the 
model point is close to the 3<r contour), while the steep spectrum 
model point lies well outside the 3<r contour. 

These results suggest that, in a statistical sense, the correla- 
tions between radio luminosity, X-ray luminosity and mass of ac- 
tive black holes require a radiatively inefficient accretion flow cou- 
pled with a (scale invariant) synchrotron emitting jet. It is important 
to stress here that this result does not suggest that the ADAF model 
for the accretion flow is the correct one, and convection and/or pow- 
erful outflows are unimportant. In fact, the nature of the radio emis- 
sion and the mere existence of the fundamental plane are clear indi- 
cations that the full dynamical model for the accretion flow should 
at some level include, in a self-consistent manner, the effects of the 
backreaction from the outflowing gas. However, our results rather 
dictates that the radiative efficiency of the flow must be low. This 
conclusion holds over the observed range of Lx/^Edd, and there- 
fore for sources which are substantially sub-Eddington. In 36.6l we 
discuss what should be expected from similar studies performed on 
samples that include a larger number of very luminous black holes. 

The jet synchrotron model for the X-ray emission with p = 2 
is only marginally consistent with the data. However, steeper elec- 
tron distributions (p = 3, for example) do predict values for the 
correlation coefficient that are well within our 3-cr contours. An as- 
sessment of the relevance of such models should be made taking in 
due account the effect of cooling, as we discuss in 36.41 



6 DISCUSSION 

The main result of our work is the discovery of a "fundamental 
plane" of black hole activity. That is, if we define the instantaneous 
state of activity of a black hole of mass M (in units of solar masses), 
by the radio luminosity at a fixed frequency (for example at 5 GHz), 
and by the hard X-ray luminosity (for example in the 2-10 keV en- 
ergy band) of its compact core, and represent such an object as a 
point in the three-dimensional space (log Lr, log Lx, log M), all 
the points representing black holes (either of stellar mass or super- 
massive) will lie preferentially on a plane, described by equation 

We note that in the case the radiatively inefficient disc model 
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(q = 2.3) the fundamental plane equation J5) implies that the radio 
luminosity satisfies: 

L R , 9=2 .3 oc rh^M 1 - 38 = M 1 ' 38 , (14) 

i.e., Lb. scales with the physical accretion rate only. This is very 
close to the predicted dependence of Lr oc til 1 ' 42 from 15.11 
for the canonical parameter choice of a = 0, p = 2, and 
<f)B oc M~ x ' 2 rh 1 ' 2 , for which the kinetic jet power Wj e t is di- 
rectly proportional to the physical accretion rate, Wj a t oc M (see 
also Falcke & Biermann 1995). Then, the total power released 
by the accretion/jet system may be written as Wtot — A/c 2 = 
Lhdl + Wj c t + V^adv.conv, where the first term on the right hand 
side is the total radiated luminosity and the last one include con- 
tributions from the energy adverted and/or stored in the convec- 
tive motions. Our results suggest that the flow must be radiatively 
inefficient, therefore, for small enough accretion rates we have 
L hol ~ Mmc 2 < tile 2 ~ Wjet + W a dv,conv On the other hand, 
Wjet oc Wadv.conv oc A/c 2 . Therefore, the issue of what the rel- 
ative fraction of the total accretion energy dissipated into the jet 
is (or, alternatively, of when a source is "jet dominated"; Fender, 
Gallo & Jonker 2003; Falcke, Kording & Markoff 2003) reduces 
to the determination of the value of the constant V^jet/Wadv.conv 
This requires the knowledge of the structure functions ip in eq. |6j, 
i.e. the specification of a jet model or the direct measure of the total 
kinetic power carried by the jet, together with a dynamical model 
for the disc-jet coupling, which is clearly beyond the scope of this 
paper. 

6.1 The thickness of the fundamental plane 

The fundamental plane defined in eq. J5J is not, however, razor thin. 
The sources are substantially scattered around it, with a dispersion 
a = 0.88 in log Lr (corresponding to a dispersion of a± = 0.62 
perpendicular to the plane). Such a scatter is not at all surprising. 
Theoretically, it can be explained in large part by the scatter in the 
radio spectral index or. For example, if we consider the ADAF 
(q — 2.3) model with fixed p = 2, we expect the theoretical rela- 
tionship Lr = (0.62+0.29a R ) log L x +(0. 80-0. 62a R ) log M+ 
&r. For the sake of simplicity, let us assume now that the unknown 
radio spectral index of all the observed sources is normally dis- 
tributed, with qr = (qr) ± cr a . In order to fit the observed 
data, we need (qr) ~ and 6r w 7. The dispersion in the ra- 
dio luminosity due to the dispersion in the unknown parameter cvr 
is then ctr > a a \/0.33 2 (A log M) 2 + 0.29 2 (A log L x /Af) 2 ~ 
2o a where 2 A log M is the range in Af and 2Alog£ x /Af the 
range in luminosity spanned by our sample. Thus, a dispersion in 
the unknown radio spectral index of « 0.3 (consistent with the ob- 
servations), can give a large contribution to the observed intrinsic 
scatter. 

In addition to the scatter produced by the diversity in o;r 
(which introduces scatter in the fundamental plane relation sim- 
ply by the fact the orientation of the plane is sliiightly different for 
different cxr) the intercept 6r of the plane also varies as we look 
at different jet models and different spectral indices. It is not clear 
a-priori how much it does so, as this depends on how the shape 
functions tp / (see section l5~l1 vary in order to produce the different 
qr. While for one spectral class with uniform qr the normaliza- 
tion cancels out when determining the correlation coefficients, this 
is not true when comparing jets from different spectral classes, i.e., 
with varying or. The fact that the scatter in the relation is moder- 
ate indicates, however, that the change in ipf over the range of radio 
spectral slopes that enter our sample is also moderate. 



Ideally, we should therefore restrict the above analysis to sub- 
classes of sources which follow a very narrow range in qr. How- 
ever, at the current stage the data do not allow such a division, both 
because the sample is too small and because the spectral indices 
are often not known to an accuracy that would allow such a treat- 
ment. The rough division of our sample into flat and steep spectrum 
sources nonetheless shows that such a treatment is possible and 
does lead to a reduction in scatter in the well defined and well sam- 
pled class of flat spectrum sources, for which we can be confident 
that the radio emission does indeed originate in the self-absorbed 
core of the jet (a prerequisite of the scale invariance model by Heinz 
& Sunyaev, 2003). 

A major source of scatter could be relativistic beaming, that 
we discuss in the following section. Also to be taken into ac- 
count is the influence of the black hole spin on jet formation 
I Bland ford & Znajek 1977] IKoide et al. 2002L While we still ex- 
pect the accretion disc to have significant influence on the jet power, 
as the magnetic field necessary to tap the black hole rotational en- 
ergy must be provided and/or confined by the accretion disc (thus 
setting the field strength and the jet power; see Livio, Ogilvie & 
Pringle 1999; Meier 2001), the spin itself will enter as a sensitive 
parameter into the radio luminosity as well, introducing scatter in 
any correlation. Once again, strong future constraints on the tight- 
ness of this relation may help constrain the possible influence of 
black hole spin (see e.g. Fender 2003, for a discussion of the GBH 
case). 

6.2 Relativistic beaming and selection effects 

Because the jets that dominate radio emission from active black 
holes are known to be relativistic, Doppler boosting of synchrotron 
radiation (relativistic beaming) must be taken into account when 
studying these sources. 

For the purpose of our study it is of primary importance to 
avoid strongly beamed sources, whose intrinsic (unbeamed) radio 
luminosity can be difficult to infer. This is the reason why BL Lac 
object have been excluded from our study. Furthermore, it is also 
important to take into account possible strong selection effects due 
to relativistic beaming, that would skew the observed correlations. 
Then, in a sample of randomly oriented jetted sources (but with- 
out the aligned ones), relativistic beaming should be an additional 
source of scatter in any correlation involving radio luminosities, the 
magnitude of which depends on the average relativistic speed of the 
jets. 

For the SMBH in our sample, the first selection criterion is the 
availability of a mass measurement of the central black hole. Thus, 
the main selection biases do not come from flux limits, and no pref- 
erential selection of bright (beamed) sources should be operating, 
although anti-selection effect might be present. 

Possible contaminations from beamed sources, however, 
could still be present among flat spectrum nearby LLAGN or type 
1 Seyfert (we have 27 of them in the sample), given the current un- 
certainties on the exact nature of their radio emission. Nevertheless, 
there are many reasons to believe that these sources are not biased 
towards jet axes close to the line of sight. Flat spectrum LLAGN 
were originally selected in the optical band, from the Palomar 
spectroscopic survey l Ho, Filippenko, & Sargent 1997a l. Many of 
those type 1 low-luminosity Seyfert and LINERS do indeed possess 
double-peaked broad emission lines l Ho et al. 2000 1, clearly indi- 
cating a non negligible inclination angle. The radio properties of a 
number of these sources were also studied by Nagar et al. (2002). 
They found that LLAGN lie on the extrapolation of the FR I/FR 
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II correlations in the Lr, — Lhc plane, while a substantial offset 
would be expected if they were all strongly affected by beaming. 
Once again, we are led to the conclusion that nearby LLAGN and 
Seyfert 1 included in our sample have a random orientation angle, 
that should disfavor strongly beamed sources. Finally we note that 
even if a few of these source were indeed strongly Doppler boosted, 
their small number will not affect significantly the measured value 
of a correlation coefficient calculated in logarithmic space (loga- 
rithmic weighting). 

Selection effects due to beaming in Galactic sources may be 
more subtle. Due to the overall scaling of radio luminosity with 
mass, GBH are expected to be on average more radio quiet than 
SMBH (see also H6.6> . It is indeed well known that galactic X- 
ray binaries are difficult to detect in the radio band, and this might 
already determine a selection bias towards beamed sources. In par- 
ticular, sources selected from GBI monitoring (see i|2.2> might be 
more beamed. If indeed GBH were all affected by strong selection 
biases toward beamed emission, the derived correlation coefficients 
of the fundamental plane, £rm in particular, might be incorrect. 
However, we know from observation of proper motion during radio 
outbursts that the jets of the classical microquasars GRO J 1655-40 
and GRS 1915+105 have large angle with the line of sight (see e.g. 
Fender 2003). Furthermore, for all the other GBH in our sample, 
apart from LS5039, orbital parameters have been inferred from op- 
tical spectroscopy, and in none of them the inferred inclination an- 
gle appears to be small. Therefore, unless the black hole jets are all 
misaligned with the orbital angular momentum I Maccarone 2002 1, 
we can exclude that strong beaming affects our sample. 

A similar conclusion had already been reached by an ac- 
curate study of the observed scatter (~ 0.7) about the Lr — 
Lx correlation in a larger sample of low/hard state GBH 
I Gallo , Fender, & Pooley 2003) . An analogous argument could be 
put forward for our entire sample, with the similar conclusion 
(given the similar amount of scatter found here) that the maximum 
Doppler factor should not be too large. This conclusion would be 
further strengthen by noting that the observed scatter about the fun- 
damental plane is not much larger than what expected from the in- 
trinsic scatter in the radio spectral index qr (see $b.\\ . However, 
precise constraints on the Doppler factor could only be obtained by 
analyzing more carefully selected samples with well determined 
(possibly fixed) radio spectral indices. 

6.3 Steep spectrum sources 

For the steep spectrum sources in our sample, the origin of the radio 
emission is not entirely clear. Steep spectra in general imply opti- 
cally thin emission. This could be an indication that the emission is 
dominated by regions far away from the central engine, e.g. in the 
diffuse extended radio lobes or hot spots. This would prohibit any 
treatment that makes use of the scale invariant model, though scal- 
ing arguments for the large scale emission still apply l Heinz J 2002t . 
Steep spectrum sources should then be excluded from the sample 
completely, and the fact that they still fit into the plane would then 
mostly be due to selection effects. It should be noted that in this 
case the radio luminosity should strongly correlate with the black 
hole mass (which determines the kinetic jet power Wj c t and thus 
the large scale radio luminosity), while it should not be related to 
the instantaneous X-ray flux, because the radio lobes reflect the 
mean power output by the central source, averaging out its tempo- 
ral variability. We expect this to be a problem for Galactic black 
hole sources, where we have lower resolving power. For supermas- 
sive black holes and for nearby LLAGN in particular, we can be 



more confident that the large scale emission is well resolved and 
thus does not contribute to the core emission used in our sample. 

It is also possible that some of the steep spectrum sources 
can be associated to core jet emission which is optically thin 
at radio frequencies, i.e., their self-absorption break lies below 
5 GHz. Because we would expect the self-absorption break fre- 
quency to be lower for higher M and lower accretion rates m 
I Heinz & Sunyaev 2003 1, this effect would imply that steep spectra 
should predominantly be observed in SMBH which arc relatively 
X-ray dim. While the trend with Lx is not clear in the data, the 
fact that the steep spectrum sources in our sample almost exclu- 
sively belong to the SMBH class does suggest that at least some 
of the sources in the sample might indeed be core dominated steep 
spectrum sources, for which the theoretical analysis of Sj5|holds. 



6.4 X-ray emission from jets 

Our analysis, taken at face value, does suggest that accretion mod- 
els for the origin of the X-ray emission fit the data somewhat better 
than pure jet models. However, at X-ray energies the effects of cool- 
ing on the particle spectrum cannot be ignored. In the context of 
scale invariant models, taking cooling into account is not possible. 
Proper treatment of particle transport including the effects of cool- 
ing may lead to different results for the synchrotron X-ray model. 
In fact, simply using p = 3 in eq. 1 1 3i (as expected in a continu- 
ous particle injection model with cooling) moves the flat spectrum 
point close to the la contour of the observed correlation in Fig. [5] 

For the case of the galactic black hole GX 339-4 in the 
hard state, where the radio-X-rays correlation is well established, 
Markoff et al. (2003) were able to fit the observed slope (£rx, 339 — 
0.71) with a synchrotron X-ray model assuming p — 2. 15 (and thus 
Qx = 0.58, which is also close to the average slope of the X-ray 
spectral index of the source in the hard state). Here we have found 
that, in a statistical sense, higher values of p are needed to explain 
the observed correlation for the entire sample in the framework of 
the jet synchrotron model for the X-ray emission. 

Furthermore, for p — 2 the X-ray synchrotron luminosity 
from the jet scales like Lx, synch oc M 1 ' 25 , while the disc X-ray 
luminosity follows Lx,disc oc M. Thus, if the X-ray emission in 
GBH jets like GX 339-4 were indeed due to synchrotron emission, 
we would expect an even stronger contribution from synchrotron 
X-rays relative to the disc emission in AGNs with the same ac- 
cretion rate, where the general wisdom (mainly based on accurate 
analysis of detailed X-ray spectra) is still that the X-rays originate 
predominantly in the disc. For p = 3, on the other hand, the syn- 
chrotron X-rays follow Lx, synch oc M, and the relative contribu- 
tion from jet and disc will be independent of M. 

Finally, the X-ray spectral variability behaviour of both GBH 
and Seyfert galaxies is consistent with pivoting of the spectrum at 
X-ray energies. If the X-ray emission is produced by synchrotron 
radiation and there is no cooling break at longer wavelengths, such 
a power-law variability would imply a radio variability which is or- 
ders of magnitude larger than observed (see discussion in Zdziarski 
et al. 2003). 

It seems, therefore, that if the X-ray emission in some of the 
sources is dominated by synchrotron emission from the jet, the ef- 
fects of cooling must play an important role in determining the ob- 
served correlation. In the context of the scale-invariant model, it 
is rather problematic to take cooling into account; this is why a 
theoretical assessment of the role of the cooling on the observed 
correlation is beyond the scope of this paper. 
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6.5 Estimating black holes masses and the nature of ULXs 

Due to the intrinsic scatter, the predictive power of the log M — 
log Lx — log Lr, relation is currently limited in scope. Nonetheless, 
from the fundamental plane equation ^5}, we can derive a formal 
relationship between the observed fluxes in the radio (at 5 GHz) 
and in the X-ray (2-10 keV) band (in erg s _1 cm -2 ), the distance 
to a source D (in Mpc) and its unknown mass. We obtain 

log M ~ 16.3 + log D + 1.28(log F R - 0.60 log F x ) ± 1.06.(15) 

In particular, we can predict in which region of log Lr — 
log Lx space we would expect to find intermediate mass black 
holes, if they exist, and specifically, where we would expect Ultra- 
Luminous X-ray sources (ULX; Makishima et al. 2000) based on 
their X-ray flux, should they be intermediate mass black holes, as 
opposed to intrinsically beamed low mass black holes. 

A clear prediction of the fundamental plane equation in this 
sense is that, if ULX are stellar mass objects (with M < 30), and 
their radio emission is not strongly Doppler boosted, there is an 
upper limit to their radio flux at 5 GHz. For example, a ULX with 
and X-ray luminosity (in the 2-10 keV band) of 10 40 ergs s -1 at a 
distance of 10 Mpc will have 

log(F R /Jy) < - 6.5 - 2 log(L>/ 10) + 0.6(log L x /40) ± 0.9(16) 

so its flux will be at most few microjansky. 

As an example, we refer to the recent claim of the observation 
of synchrotron radio emission from the ULX 2E 1400.2-4108, in 
the dwarf irregular galaxy NGC 5408 (at a distance of ~ 4.8 Mpc) 
jKaaret et al. 2003t . The observed radio (0.26 mjy at 4.8 GHz) and 
X-ray fluxes (2.7 x 10" 12 erg s' 1 cm' 2 in the 0.3-8 keV band), 
if indeed both associated with a compact source at the distance of 
NGC 5408, would imply a mass of log M ~ 4 ± 1. However, 
we should stress here that such an estimate only holds the radio 
emission is not beamed. The uncomfortably high value for the esti- 
mated mass, the steep inferred radio spectral index and the radio to 
X-ray flux ratio, all suggest that the observed intensity of the radio 
emission in this source may be Doppler boosted by a relativistic 
jet pointing in our direction as already suggested by Kaaret et al. 
(2003). 

6.6 Accretion mode changes 

It is well accepted, both from theory and observations, that ac- 
cretion can proceed in different modes (or states), with different 
radiative efficiencies and spectral properties l Abra mowicz T9981 
IDone 20021 prank, King, & Raine 2002) |Narayan 200*2} , mainly 
driven by variations of the dimensionless accretion rate rh. Here 
we have demonstrated (SJ5J that the correlation between radio and 
X-ray luminosity in GBHs and SMBHs can provide valuable con- 
straints on the emission mechanisms in these bands and on the 
physics of accretion. We have shown how low radiative efficiency 
accretion models can possibly provide the underlying scalings that 
most easily explain the properties of the observed fundamental 
plane. Radiatively efficient thin disc-corona systems are clearly 
inconsistent with the data. The only possibility for these mod- 
els to be reconciled with our observations would entail a highly 
radiatively inefficient corona, where most of the dissipated mag- 
netic energy is converted into bulk kinetic energy of outflowing 
gas l Beloborodov 1999 1, as described by Merloni & Fabian (2002). 
This would require some specific model for the magnetic dissipa- 
tion processes, and for their scaling with M and m, which are cur- 
rently unknown. 



However, because both such modes of accretion are expected 
to occur only below accretion rates < few percent of Edding- 
ton, we would expect the log M — log Lr — log Lx correlation 
to break down at around this limit. In other words, we should ex- 
pect the QSOs and the bright Seyferts in our sample, which occupy 
the region of high accretion rates and are independently known to 
have spectral characteristics inconsistent with models of low ra- 
diative efficiency, to depart from the observed correlations. Also 
in GBHs, it has indeed been shown that the correlation between 
radio and X-ray luminosity breaks down as the sources switch to 
their high states <Maccarone 2003l[Gairo, Fender, & Pooley 200'3') . 
In particular, high luminosity states of GBH (high or very high 
states) show highly variable radio fluxes and radio to X-ray ratios 
(Mirabel & Rodr iguez 1994| IFender et al. 19991 . If this reflects a 
general property of the disc-jet coupling at high accretion rates, a 
substantial increase in the scatter should appear above the critical 
accretion rate. 

In order to illustrate this point further, in FigureQwe plot, as 
a function of the ratio Lx/L-^dd, the radio luminosity divided by 
M 1 ' 38 (such a scaling for the radio luminosity with mass is ob- 
tained directly from the fundamental plane equation {5J by impos- 
ing that the X-ray luminosity scales linearly with black hole mass). 
As indeed expected, by rescaling the radio luminosity in such a 
way all the different tracks corresponding to different mass bins in 
Fig.|3j) collapse into a single one (with some residual scatter). The 
region between the two vertical lines corresponds to the expected 
values of Lx/L^dd above which a change of accretion mode, from 
radiatively inefficient to standard radiatively efficient is expected to 
occur. Our SMBH sample is still limited in order to test whether 
such a change is observed. In fact, we select out of our sample the 
majority of bright quasars, both radio quiet and radio loud, because 
of the lack of a reliable mass estimates for these sources. However, 
there is indeed a hint of an apparent increase in the scatter about the 
correlation, above the critical accretion rates, as seen in Figure[7| 

We may thus speculate that the famous (and still much de- 
bated, see e.g. Cirasuolo et al. 2003) radio loud/radio quiet di- 
chotomy of quasars will appear only at the highest values of m, 
and be caused mainly by a switch of accretion mode analogous to 
the high/very high transition in GBH. At low accretion rates, black 
holes seem to follow the more regular behavior circumscribed by 
the fundamental plane of eq. J5J. Such sources not only tend to 
be radio loud jHo & Peng 200 1| IHo 2002at , but also their radio 
loudness parameter, Fx (here defined as the ratio of radio to X- 
ray luminosity), obeys the following scaling: Rx = Lr,/Lx oc 
L^ 2//5 M 4//5 . Therefore, the smallest the X-ray luminosity, the 
more radio loud these sources are tHo 200 2a I. In this regime, no 
dichotomy has to be expected, as already suggested by Nagar et al. 
(2002). 



7 SUMMARY 

We have shown how the coupling between accretion flows and 
jet emission in black holes can be understood by studying large 
samples of both galactic (GBH) and supermassive (SMBH) black 
holes with measured mass, observed at both radio and X-ray wave- 
lengths. 

We have compiled, from the existing literature, a sample of 
~100 supermassive black holes and ~50 different observations of 
galactic, stellar mass black holes. We required that the SMBH have 
a reliable measure of the central black hole mass, either (preferen- 
tially) direct, via high resolution kinematics studies of surrounding 
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Figure 7. The radio luminosity logLj^, divided by M 1,38 as a function 
of the ratio Lx/^Edd- Solid squares are for steep spectrum sources, solid 
triangles for flat spectrum ones and open circles for sources with undeter- 
mined radio spectral index. Two vertical lines mark the boundary of the re- 
gion where we expect the critical luminosity for the mode change between 
radiatively inefficient and efficient accretion. The scaling for the radio lumi- 
nosity with mass is obtained directly from the fundamental plane equation 
l5l by imposing that the X-ray luminosity scales linearly with black hole 



possible to predict the correct amount of scatter for any such rela- 
tionship. 

By comparing the observationally derived correlation coeffi- 
cients to the theoretically predicted ones, we are able to put con- 
straints on accretion models and on the disc-jet coupling. We 
demonstrate that the X-ray emission from black holes accreting at 
less than a few per cent of the Eddington rate cannot be produced 
by radiatively efficient accretion, while radiatively inefficient ac- 
cretion flows agree well with the data. The optically thin X-ray jet 
synchrotron emission model is only marginally consistent with the 
observed correlation; however, cooling of the electrons should be 
properly taken into account before drawing a firmer conclusion on 
the issue. 

The fundamental plane equation also implies that low lumi- 
nosity black holes should be more and more radio loud as the ac- 
cretion rate decreases, as observed I Ho 2002a I. The possible di- 
chotomy between radio quiet and radio loud Quasars could in- 
stead be due to a switch of accretion mode at the highest accretion 
rates. Therefore, more accurately selected samples of more lumi- 
nous sources, including a significant number of both radio quiet 
and loud quasars with well determined radio spectral properties, 
should allow the observation of bends, discontinuities or even bi- 
furcations of the fundamental plane associated with accretion mode 
switches and with the turning on and off of the most powerful radio 
sources. 



stars and gas or reverberation mapping of the broad line region (for 
luminous Seyfert 1 galaxies and Quasars), or, indirectly, thanks to 
the empirical correlation between black hole masses and central ve- 
locity dispersion. Although such selection criteria introduce a num- 
ber of biases in the final distribution of M, which are difficult to 
control, our final sample spans a large enough range in both mass 
and accretion rate to unveil major trends in the physical properties 
of active black holes. This is a crucial point, and represents the main 
novelty of our approach. Most previous studies on the subject have 
only considered correlations of jet (radio) luminosity with black 
hole mass or with luminosity in other wavewbands (and hence with 
accretion rate) separately. 

We have performed a partial correlation analysis on the sam- 
ple, and concluded that the radio luminosity is strongly correlated 
with both black hole mass and X-ray luminosity, while in turn, the 
X-ray luminosity correlates with both mass and radio luminosity. 
In light of this, we have carried out a multivariate linear regression 
analysis and derived the best fit correlation coefficient between ra- 
dio luminosity, X-ray luminosity and black hole mass. 

The data from the entire sample allow us to put tight 
constraints on the correlation coefficients. Remarkably, we find 
that the sources lie preferentially on a plane ("the funda- 
mental plane of black hole activity") in the three dimen- 
sional (log Lb,, log Lx, log M) space, described by the equation 
log Lr = (0.60t°;ii)logL x + (0.78™) log M + 7.33t£$. 
The scatter around this plane is, however, significant (<tr = 0.88). 

From a theoretical point of view, these results clearly sug- 
gest that the ansatz of scale invariance for the disc-jet coupling 
jHeinz &"Su nyaev 2003 1 captures the main physical properties of 
such systems. Thus, a universal theoretical scaling between the ra- 
dio flux at a given frequency and both mass and accretion rate can 
be derived, independently of the jet model, with scaling indices that 
depend only on the (observable) spectral slope of the synchrotron 
emission in the radio band, and on the accretion mode. Also, it is 



NOTE ADDED IN PROOFS 

After the submission of this manuscript another work has appeared 
as preprint I Falcke, Kording, & Markoff 2003 l, that is similar in 
scope, and reaches very similar conclusions about the scale invari- 
ant properties of low luminosity black holes. 
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